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1  Technical  Objectives 

Coanda-effect  circulation  control  (CC)  lias  shown  tremendous  promise  for  liigli-lift  control  sur¬ 
faces  in  underwater  naval  vehicles  (Rogers  &  Donnelly  2004).  CC  increases  lift  generation  at 
low  mean  flow  speeds  without  deploying  mechanical  devices.  However,  additional  self-noise  is 
produced  over  and  above  trailing  edge  noise.  Unfortunately,  our  understanding  of  CC  acoustics 
is  limited.  Salikuddin  ei  al.  (1987)  and,  more  recently,  Munro  et  ai  (2001)  have  experimentally 
studied  the  noise  due  to  conventional  upper  surface  blowing  on  an  airfoil  with  a  mechanical 
flap  appendage,  which  is  most  suitable  for  aircraft  configurations.  On  the  other  hand,  Howe 
(2002)  has  analyzed  sound  production  from  a  hydrofoil  and  identified  three  uiecliaiiisms:  (1) 
low  frequency  curvature  noise  associated  with  interaction  of  a  turbulent  boundary  layer  with 
the  rounded  trailing  edge,  (2)  high  frequency  jet-slot  noise,  and  (3)  passive  slot  noise  caused 
by  diffraction  of  near  field  pressure  fluctuations  in  the  turbulent  boundary  layer  by  the  slot  lip. 
This  theory  has  not  been  experimentally  confirmed.  To  complicate  matters  further,  recent  w'ork 
by  Rogers  t  Donnelly  (2004)  has  demonstrated  the  potential  of  dual-slotted  (i.e..  upper  and 
lower  surface)  blowing  configurations.  The  noise  associated  with  this  configuration  has  not  yet 
been  systematically  studied. 

The  objectives  of  this  study  are  to  experimentally  characterize  the  acoustics  and  fluid  dy¬ 
namics  of  an  elliptic  circulation  control  airfoil.  Particle  image  velocimetry  measurements  provide 
insight  into  the  flow  physics  related  to  the  noise  generation  mechanisms.  A  variety  of  acoustic 
measurement  techniques  are  compared  to  determine  what  methods  may  be  suitable  for  mea¬ 
suring  circulation  control  noise.  Noise  sources  are  identified,  and  their  far-field  acoustic  levels 
are  compared  to  determine  the  primary  sources  over  different  frequency  ranges.  Finally,  Hewe's 
theoretical  model  of  circulation  control  noise  is  evaluated  and  compared  with  measurements 
(Howe  2002). 

2  Technical  Approach 

A  two-dimensional,  dual-slotted,  elliptic  circulation  control  airfoil  based  on  the  hydrofoil  studied 
by  Rogers  Sz  Donnelly  (2004)  is  designed  and  fabricated.  Tests  are  conducted  in  the  University 
of  Florida  Aeroacoustic  Flow  Facility  (UFAFF)  at  low  and  moderate  chord  Reynolds  numbers 
{Rcr  =  6.5  X  10^  to  1.3  X  10®)  and  zero  angle  of  attcvck.  Focus  is  placed  on  low  blowing,  moderate 
lift  performance  typical  of  an  underwater  vehicle  control  surface.  Higher  lift  cases  appropriate  for 
application  to  STOL  aircraft  are  also  studied.  Steady  surface  pressure  measurements  are  us(‘(i  to 
characterize  the  freestream  flow  and  determine  lift  as  a  function  of  moiueiituiu  (‘oeflicient.  The 
freestream  and  curved  wall  jet  flows  are  measured  and  analyzed  using  PIV.  Multiple  far- field 
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microphones  and  a  nested  phased  microphone  array  are  used  for  determining  source  levels  and 
identifying  noise  sources.  Using  data  from  the  PIV  study,  Howe's  model  of  circulation  control 
acoustics  is  evaluated  and  compared  with  measurements  (Howe  2002). 

3  Experimental  Setup 

For  a  complete  description  of  the  experimental  setup,  please  refer  to  Wetzel  (2011).  For  conve¬ 
nience,  a  brief  summary  of  the  experimental  setup  is  included  in  this  report. 

3.1  Circulation  Control  Airfoil 

A  dual-slotted,  two-dimensional  circulation  control  airfoil,  shown  in  Figure  lA,  is  designed  based 
on  the  geometry  of  the  hydrofoil  previously  studied  by  Rogers  &  Donnelly  (2004).  The  airfoifs 
profile  is  a  20%  thickness-to-chord  ratio  ellipse  with  a  cylindrical  trailing  edge  and  no  camber. 
Details  regarding  the  trailing  edge  geometry  are  providwl  in  Figure  IB.  The  airfoil  has  a  0.5207 
m  chord  and  1.12  m  span.  Although  the  frontal  area  ratio  between  the  model  and  tunnel  inlet 
cross-sections  is  somewhat  high  (14%),  the  large  airfoil  size  is  chosen  for  a  number  of  reasons, 
including:  to  suit  anticipated  full-scale  application  chord  Reynolds  numbers  and  inomentum 
coeflicieiits,  for  improved  spatial  resolution  in  measureinents,  for  reduced  coinpressibilily  eli'ecls 
in  the  slot  jet,  and  for  lower  acoustic  frequencies.  The  height  of  each  blowing  slot  is  adjustable 
using  eight  sets  of  push-pull  screws  equally  spaced  across  the  span.  Each  blowing  slot  is  supplied 
pressurized  air  from  two  constant  area  air  feeds,  one  on  each  side  of  each  plenum.  The  model 
is  primarily  machined  from  cast  aluminum  tool  plate. 

During  the  investigation,  only  a  single  blowing  slot  is  used  at  any  given  time.  The  other, 
unused  blowing  slot  is  set  to  the  nominal  slot  height  of  1.0  mm  and  sealed  with  tape.  Spaiiwise 
strips  of  12  mm  wide,  0.4  nun  thick,  Glasfaser  Flugzeug-Service  zig-zag  turbulator  tape  are 
placed  on  the  airfoil  at  18%  chord  to  trip  the  upper  and  lower  surface  boundary  layers,  ghboring 
labs  for  model  leak  tests  and  slot  flow  uniformity  studies  between  tunnel  entries. 

3.2  University  of  Florida  Aeroacoustic  Flow  Facility 

The  UFAFF  is  an  open-return  wind  tunnel  with  an  open-jet  test  section  installed  in  an  ISO 
3745-cert ified  100  Hz  anechoic  chaml)er,  which  is  used  to  simulate  a  free  field,  or  reflection- 
free  acoustic  environment.  The  test  section  measures  0.74  m  tall  by  1.12  m  wide  and  extends 
1.83  m  in  the  flow  direction.  Test  section  speeds  between  18  and  75  m/s  can  be  reached  with 
turbulent  intensity  levels  below  0.1%  at  frequencies  above  10  Hz  (Mathew  et  al.  2005).  hi  2010. 
a  flow  silencer  was  added  upstream  of  the  inlet,  and  over  one  meter  of  the  diffuser  was  rcunoved. 
Figure  2  provides  schematics  of  the  circulation  control  airfoil  installed  in  the  pre-  and  post¬ 
modification  tunnel  configurations.  For  illustration  purposes,  some  ceiling  wedges  are  removed 
in  the  renderings  to  reveal  the  chamber  interior  and  test  section.  The  airfoil  is  installed  13 
cm  downstream  of  the  inlet  exit  plane  at  zero  degrees  geometric  angle  of  attack.  The  model 
is  always  bounded  by  sidewalls  for  two  diinensional  flow.  For  more  iiifornmtion  regarding  llu' 
UFAFf\  please  refer  to  Mathew  et  al.  (2005). 
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3.3  Flow  Measurements 


3.3,1  Steady  Pressure  and  Lift 

Steady  surface  pressure,  is  measured  using  47  0.711  mm  inner  diameter  midspan  static  pres¬ 
sure  ports,  including  nine  located  on  the  midspan  trailing  edge  instrument  ring.  The  taps  are 
symmetrically  distributed  on  the  upper  and  lower  surfaces  of  the  airfoil.  Two  taps  place<l  at 
the  leading  edge  of  the  airfoil  are  slightly  offset  from  midspan  due  to  installation  constraints. 
Three  16-channel  Esterline  Pressure  Systems  (PSl)  pressure  scanners  (ranges  5  psi,  1  psi,  and 
10  inH20)  measure  39  static  surface  pressures,  four  Venturi  meter  pressures,  and  the  blowing 
slot  plenum  pressure.  Data  are  acquired  at  two  samples  per  second  for  30  seconds.  The  pres¬ 
sure  scanners’  reference  pressure  ports  are  connected  to  the  tunnel  pitot-static  pressure  port. 
The  pitot-static  and  pitot-dynamic  pressures  are  measured  using  Mensor  6115  digital  pressure 
transducers  with  ranges  of  16  psi  a  and  0.58  psig,  respectively.  Surface  pressure  coefficients  are 
computed  by  dividing  the  pressure  scanner  output,  Ps  l>y  the  pitot-dynamic  pres.sure, 


c„  = 


Ps  -  P  x 
Qoc 


(1) 


Because  of  the  coarse  resolution  of  pressure  taps,  the  Cp  data,  plotted  versus  normalized  chord- 
wise  position,  are  fit  with  a  cubic  spline  curve  that  is  integrated  to  determine  the  lift  coefficient 
per  unit  span. 

Cl  —  I  ^Cpd{x/c)  =  /  {(^p.lower  ”  ^ p,uppf"r)  d{x j c)  (2) 

Jo  Jo 


3.3.2  Slot  Jet  Velocity  and  Momentum  Coefficient 


The  mean  slot  jet  velocity  is  computed  assuming  an  isentropic  expansion  from  the  plenum  to 
the  frees tream, 


Ujet  = 


;>o  / 


(7-l)/7" 


(3) 


where  7  is  the  ratio  of  specific  heats,  and  R  is  the  ideal  gas  constant.  The  measured  fn‘estream 
pitot-static  pressure,  is  added  to  the  measured  relative  plenum  pressure  to  determine  the 
absolute  plenum  pressure,  po.  The  air  line  stagnation  temperature,  Tq,  is  recorded  using  a 
resistance  temperature  detector  (RTD).  For  the  complete  derivation  of  Equation  3  and  the 
influence  of  a  nonzero  plenimi  velocity  on  the  je't  velocity  estimate,  see  Wetzel  (2011). 

The  momentum  coefficient  is  calculated  using  the  combined  mass  flow  rates  of  the  two 
Venturi  meters,  the  jet  velocity  computed  from  Equation  3,  the  measured  freestream  dynamic 
pressure,  and  the  circulation  control  model  planform  area.  5, 


ihUjet 
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3.3.3  Particle  Image  Velocimetry 

Particle  image  velocimetry  provides  a  non-intrusive  means  of  measuring  the  velocity  and  relatecl 
properties  of  a  flow  field  (Raffel  et  ai  2007).  The  local  displacement  of  tiny  particles  seeded 
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within  a  flow  are  detecttHi  by  imaging  two  snapshots  of  the  flow  ilhiniiiiatod  by  a  “instantaneous” 
light  source.  The  light  source,  which  is  a  laser  beam  that  passes  through  a  series  of  lenses  to  form 
a  thill  laser  sheet,  illuniiiiates  a  planar  region  of  interest  in  the  flow  field.  The  two  spatially 
coincident  lasers  are  pulsed  at  a  known  rate  but  separated  by  a  specified  amount  of  time', 
typically  on  the  order  of  microseconds.  For  two-dimensional  PIV,  a  camera  positioiuxi  outside 
the  flow  field  and  aligned  normal  to  the  laser  sheet  captures  a  pair  of  images  for  each  pair  of  laser 
pulses.  The  image  pair  is  discretized  into  smaller  regions  called  interrogation  windows,  where 
each  interrogation  window  is  designed  to  contain  several  common  illuiniiiated  particles  for  each 
frame  of  the  image  pair.  The  spatial  cross-correlation  of  eaeh  corresponding  pair  of  windows  is 
then  coinpiited  using  Fast  Fourier  Transforms  (FFTs).  The  average  x-  and  y-displaceinents  of 
the  group  of  particles  within  each  window  pair  are  determined  from  the  location  of  the  highest 
correlation  peak.  The  in-plane  velocity  components  eire  calculated  from  these  displacements 
and  the  known  time  between  laser  pulses.  Altogether,  the  discretized  velocity  ineasureinents 
obtained  from  an  image  pair  provide  a  snapshot  of  the  velocity  flow  field.  With  a  large  enough  .set 
of  measured  vector  fields  such  that  statistical  convergence  is  achieved,  turbulent  flow  quantili(?s 
can  also  be  estimated. 

Five  different  PIV  experiments  are  performed  to  capture  different  flow  pheiionienon.  The 
flow  is  ineasiircd  in  regions  around  the  leading  edge,  behind  the  trailing  edge,  above  the  slot 
lip,  and  along  the  curved  trailing  edge  in  the  chordwise  plane,  as  shown  in  Figure  3.  The 
trailing  edge  flow  is  also  measured  once  in  the  spanwise  plane.  For  more  information  on  the 
PIV  experimental  setup,  please  see  Wetzel  (2011). 

3.4  Acoustic  Measurements 

Acoustic  measurement  techniques  utilized  in  this  study  include  free-standing  microphones  and 
array-based  methods.  Both  are  discussed  in  this  section. 

3.4.1  Free-Standing  Microphones 

Two  different  sets  of  free-standing  microphone  data  are  acciuired.  The  first  is  measured  prior 
to  the  shortening  of  the  UFAFF  diffuser  and  the  addition  of  an  inlet  silencer  but  is  included 
in  this  analysis  because  multiple  slot  heights  are  tested.  For  that  experiment,  four  CI.R.A.S. 
Sound  and  Vibration  (G.R.A.S.)  6.35  mm  40BE  free-field  microphones  are  placed  outside  the 
wind  tunnel  open  jet  shear  layers.  The  microphones  are  oriented  in  pairs,  tw'o  above  and  two 
below  the  airfoil.  a.s  shown  in  Figure  4.  The  G.R.A.S.  microphones  are  calibrated  using  a  Briiel 
iz  KJa'r  (B&K)  4231  Sound  Calibrator  at  a  frecpieiicy  of  1  kHz.  G.R.A.S.  40BE  microphones  are 
designed  for  a  flat  frequency  response  for  0^  incidence  with  the  protection  grid  installed.  The 
frequency  responses  of  all  G.R.A.S.  microphones  used  deviate  by  less  than  1  dB  (ref.  20// Pa), 
so  no  corrections  are  applied  to  the  measured  spectra. 

The  second  set  of  data  are  recorded  simultaneously  with  the  phased  array,  which  is  described 
in  Section  3.4.2.  Three  microphones  are  placed  above  the  trailing  edge,  as  illustrated  in  Figure 
5.  The  center  microphone,  labeled  “'^12"  in  Figure  5,  is  a  G.R.A.S.  40BE,  while  the  upstream 
(“Ml")  and  downstream  (“MS”)  microphones  displaced  10.2  cm  from  the  center  inicrophoiie 
are  B&K  4939  6.35  mm  free-field  microphones.  All  three  microphones  are  approximately  1.12 
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m  from  the  trailing  edge.  The  microphones  arc  calibrated  using  a  B&K  4228  Pistonphone  at 
a  frequency  of  250  Hz.  The  microphones  are  designed  for  a  flat  frequency  response'  for 

0°  incidence  with  the  protection  grid  removed.  Since  the  microphone  grids  remain  installed, 
the  frequency  response  of  each  B&K  microphone  is  used  to  adjust  the  measured  spectrum 
accordingly. 

Various  processing  techniques  are  applied  to  the  microphone  data,  including  single  mi¬ 
crophone  autospectra,  coherent  output  power,  and  the  three-microphone  method  (Brooks 
Hodgson  1981;  Hutcheson  Brooks  2002;  Bahr  et  ai  2008:  Bendat  &  Piersol  2000).  These 
techniques  are  described  in  detail  by  Wetzel  (2011). 

3.4.2  Phased  Array 

Phased  microphone  arrays  are  used  extensively  for  noise  source  as.sessment  in  aeroacoustic  aj)- 
plications.  A  jdiased  microphone  array  consists  of  a  significant  number  of  niicrophoiK‘s  strato 
gically  located  in  a  frame  or  flush- mounted  in  a  flat  plate.  Tlu'  array  is  “focused”  to  a  point  of 
interest  by  appropriately  weighting  and  delaying  the  microphones’  signals  in  a  process  known 
as  beamforming. 

Beamforming  The  simplest  and  most  straightforward  beamformer  is  the  delay-and-sum 
beainfonner.  Consider  the  illustration  in  Figure  6  of  an  array  comprised  of  m  =  1,...,;\/ 
microphones  and  a  scanning  region  of  /  =  1,...,L  points  located  at  x/,  each  assumed  to  be  an 
incoherent  and  statistically  independent  inonopole  source.  The  estimated  acoustic  power  at 
each  location  is  given  by  the  following  equation  5  (Dougherty  2002;  Yardibi  ct  al.  2010a). 

Pi  =  -^af'Gai  (5) 

The  weight,  or  steering,  vector  a/  compensates  for  the  different  propagation  distances  traveled 
by  the  sound  wave  from  the  scanning  point  xi  to  the  different  microphones  and  is  defined 
by  Equation  6.  The  superscript  in  Equation  5  corresponds  to  the  Herniitian,  or  complex- 
conjugate  transpose  of  this  matrix. 
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a/  =  — 

n.o 


ri  i€ 


(6) 


In  Equation  6,  the  distance  from  each  scanning  point  /  to  each  microphone  in  is  designated 
and  the  distance  from  the  array  center  to  the  Ith  scanning  point  is  denoted  r/  Q.  Note  that 
the  1/7*/, 0  factor  is  used  to  scale  the  power  measured  by  each  microphone  to  match  the  power 
that  would  be  measured  at  the  array  center.  Finally,  G  represents  the  estimated  cross- spectral 
matrix  (CSM)  of  the  measured  microphone  signals. 

Integrated  spectra  can  be  obtained  by  summing  the  powers  of  each  scanning  point  and 
normalizing  by  the  point-spread-function  (PSF)  of  the  array  (Oelennans  ct  al.  2007;  Yardibi 


et  al.  2010b). 
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The  PSF  is  the  theoretical  array  response  to  a  inonopole  placed  above  the  center  of  the  array 
at  the  same  height  as  the  scanning  plane.  In  computing  the  power  via  Eejuation  7,  only  selected 
scanning  points  in  a  particular  region  of  interest  need  to  be  evaluated,  providing  a  means  to 
‘‘remove’’  noise  sources  other  than  the  source(s)  of  interest. 

Compared  to  the  previously  described  methods,  phased  ai'ray  measurements  arc  iniich  more 
expensive  from  cost,  time,  and  computational  standpoints.  In  addition,  while  relative  sound 
levels  from  array  integration  can  typically  be  trusted,  absolute  levels  iiiay  be  erroneously  low 
due  to  coherence  loss  between  microphones  in  the  array,  particularly  at  higher  frequencies 
(Oelermans  et  al.  2007).  However,  unlike  the  conventional  microphone  techniques  discussed 
earlier,  phased  arrays  are  extremely  useful  in  identifying  noise  sources  in  a  complex  acoustic 
field. 


Array  setup  A  60-clement  free-ficld  nested  microphone  array  is  designtxl  and  built  for  the 
circulation  control  experiments.  The  inner,  7.62  cm  aperture,  25-'element  array  is  designed  for 
frequencies  between  10  kHz  and  80  kHz  and  is  comprised  of  10  B&K  4954 B  and  15  G.R.A.S. 
40BE  high-frequency  microphones.  The  outer,  73.7  cm  aperture,  40-elcmcnt  array  is  desigiuxi 
for  frequencies  below  10  kHz  and  is  comprised  of  35  B&K  4958  “array”  microphones  in  addition 
to  five  shared  G.R.A.S.  40BE  microphones  also  part  of  the  inner  array. 

The  microphones  are  installed  on  a  stereolithography  frame  with  built-in  microphone  hold¬ 
ers.  The  frame  is  wrapped  in  Nomex  and  foam  to  mitigate  reflections  caused  by  the  frame  and 
its  mount.  The  array  is  installed  in  the  anechoic  wind  tunnel  1.12  m  below  the  airfoil’s  trailing 
edge  and  centered  at  niidspan,  as  shown  in  Figure  5.  A  picture  of  the  test  setup  along  with 
pictures  of  the  array  are  provided  in  Figure  7.  Design  microphone  coordinates  and  the  serial 
numbers  of  all  microphones  used  are  listed  in  Wetzel  (2011). 

Calibration  The  goal  of  array  calibration  is  to  characterize  the  magnitude  and  phase  rt^poiise 
of  each  microphone  relative  to  the  array  center,  or  reference,  microphone.  The  frequency  re¬ 
sponse  between  each  sensor  and  the  reference  is  diffeTOiit  due  to  the  individual  inicrophoiu' 
responses  and  scattering  from  neighboring  microphones  and  the  sensor’s  own  protective  grid. 
The  calibration  procedure  used  in  this  investigation  is  described  by  Wetzel  (2011). 

3.4.3  Data  Acquisition 

All  acoustic  data  are  acquired  at  201,800  Hz  for  30  seconds.  During  processing,  the  data  from 
each  channel  are  separated  into  480  blocks  of  12,800  samples  each,  resulting  in  a  16  Hz  bin 
width.  The  data  are  then  proccssctl  using  a  Hanning  window  and  75%  overlap,  producing  996 
effective  averages  and  a  normalized  autospectral  random  uncertainty  of  3.2%. 
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4  Fluid  Dynamics 


Even  though  the  primary  motivation  for  this  research  is  acoustics,  the  fluid  dynamics  of  the 
problem  are  equally  important.  Knowledge  of  the  flow  behavior  is  crucial  to  understaiidiiig  the 
nature  of  the  flow  and  tlicxjretical  uoise  source's,  lik('  those  outlined  by  Howe  (2002).  Surface 
prc'ssurc  ineasuremeiits  are  used  to  provide  insight  into  the  flow  around  the  airfoil.  PIV  is 
used  to  measure  details  of  the  flow  field,  including  leading  edge  stagnation  point  movement  and 
streamline  curvature,  the  influence  of  the  blowing  slot  on  the  trailing  edge  wake,  characteristics 
of  the  turbulent  boundary  layer  passing  over  the  slot  lip,  and  curved  wall  jet  development 
and  separation.  Before  the  flow  field  is  studied  extensively,  the  circulation  control  airfoil  is 
characterized. 

4.1  Freestream  Flow  Characterization 

The  freestream  flow  field  is  characterized  using  a  combination  of  static  surface  pressure  and 
PIV  measurements.  The  results  from  these  tests  are  presented  in  this  section. 

4.1.1  Surface  Pressure  Measurements 

Static  surface  pressure  distributions  are  plotted  in  Figure  8  for  throe  different  inomeiUum  (oef- 
ficients,  =  0,  0.015,  and  0.057  at  Rcc  =  6.5  x  10^  and  h/c  =  0.0019.  For  clarity,  uncertainty 
bounds  are  not  included,  as  they  are  tyj)ically  smaller  than  the  datum  point  size.  The  effect 
of  blowing  on  the  Cp  distribution  is  apparent  by  the  large  trailing  edge  suction  peak  whose 
magnitude  increases  substantially  with  larger  values  of  Cp.  A  leading  edge  suction  peak  typ¬ 
ically  observed  in  cireulation  eontrol  airfoil  experiments,  including  those  of  Abramson  (1975). 
is  absent.  Furthermore,  the  lift  produced  by  the  airfoil  is  substantially  lower  than  the  lift  pro¬ 
duced  by  Abramson’s  similar  elliptic  circulation  control  airfoil.  Lift  coefficients  for  a  variety  of 
momentum  coefficients  are  compared  in  Figure  9.  It  is  likely  that  the  absence  of  the  leading 
edge  suction  peak  is  the  source  of  the  lift  deficiency  observed  in  Figure  9. 

4.1.2  Test  Section  Influence 

To  determine  if  the  airfoil  location  in  the  test  section  is  responsible  for  the  uiissing  leading  edge 
suction  peak,  the  Jiirfoil  position  is  varied  by  13  cm  vertically  (0.25c)  and  38  cm  (0.73(')  in  the 
streainwise  direction.  However,  the  Cp  distributions  remain  unchanged.  Since  the  tests  per¬ 
formed  by  Abramson  (1975)  took  place  in  a  closed  wind  tunnel,  the  absence  of  the  leading  edge 
suction  peak  could  be  due  to  the  UFAFF  open  jet  test  section.  To  investigate  this  phenomenon 
further,  the  wind  tunnel  test  section  is  enclosed  with  foam  walls,  as  shown  in  Figure  10. 

Cp  data  measured  in  the  closed  configuration  are  corrected  to  account  for  solid  blockage 
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(Pope  &:  Rae  1984).  Specifically,  Cp,  Cp,  and  Rcc  arc  corrected  using 


c  ~  ^ 

(1  +  2.,) 

(«) 

c  - 

(1+2.,) 

(9) 

Rec,c  =  -h  f^) 

(10) 

where  tt  —  0.0357  is  one-quarter  of  the  ratio  between  the  model  frontal  area  and  test  section 
cross-sectional  area.  Once  again,  for  clarity,  uncertainty  bounds  are  not  included,  as  they  are 
typically  smaller  than  the  datum  point  size.  Figure  11  compares  niidspaii  Cp  distributions 
measured  in  both  open  and  closed  test  sections  for  Cp  =  0.  Both  Cp  distributions  have  a 
similar  shape,  but  the  closed  test  section  Cp  distribution  is  noticeably  offset  from  its  open  test 
section  counterpart.  Notably,  the  closed  test  section  Cp  data  are  more  negative,  indicative  of 
the  increased  flow  speed  caused  by  the  ciddition  of  a  floor  and  ceiling  and  their  containment  of 
the  tunnel  jet. 

Figure  12  compares  open  and  closed  test  section  Cp  distributions  for  Cp  =  0.057.  Both  Cp 
distributions  have  similar  prominent  trailing  edge  suction  peaks  associated  with  the  trailing  edge 
jet.  However,  at  the  leading  edge,  the  Cp  distributions  arc  drastically  different.  Tlu'  loading 
edge  suction  peak  absent  in  the  open  t(^st  section  configuration  appears  when  the  lest  section 
is  fully  enclosed.  Lift  coefficients  for  the  closed  test  section  and  open  test  section  cases  are 
compared  with  data  from  Abramson  (1975)  in  Figure  13.  The  closed  test  section  lift  curve  is  in 
close  agreement  with  Abramson’s  lift  values  (Abramson  applied  the  same  blockage  corrections 
used  ill  the  present  analysis).  Therefore,  the  wind  tunnel  open  jet  test  section  is  responsible  for 
the  disappearance  of  the  leading  edge  suction  peak. 

Since  the  primary  motivation  for  this  investigation  is  acoustics,  it  is  concerning  that  the 
Cp  distributions  measured  with  an  open  and  closed  test  section  differ  significantly.  The  mid- 
frequency  circulation  control  noise  source  theorized  by  Howe  (2002),  passive-slot  noise,  is  pro¬ 
duced  by  fre('stream  boundary  layer  turbulence  scattering  off  the  slot  lip.  If  the  charac  teristics 
of  the  turbulent  boundary  layer  near  the  slot  are  highly-dependeiit  on  the  flow  upstream  of  the 
slot,  then  so  too  is  the  sound  produced.  PIV  is  used  to  further  investigate  the  influence  of  the 
test  section  on  the  flow  field,  including  the  freestream  boundary  layer  passing  over  the  lip. 

4.1.3  Closed  Test  Section  Behavior 

The  leading  edge  flow  field  and  the  boundary  layer  flow  passing  over  the  slot  lip  are  measured 
using  PIV.  Please  refer  to  Wetzel  (2011)  for  information  regarding  the  experimental  sc'tup.  In 
order  to  image  the  flow  field,  one  foam  sidew^all  is  replaced  with  a  clear  polycarbonate  panel. 
The  PIV  cameras  are  mounted  behind  this  panel  outside  the  flow.  When  the  test  sc'ction  is 
fully  enclosed,  the  foam  ceiling  is  also  replaced  with  a  clear  polycarbonate  panel  to  permit  the 
laser  sheet  to  pass  through  and  illuminate  the  flow.  However,  the  tunnel  floor  remains  foam.  A 
picture  of  the  test  section  enclosed  for  PIV  measurements  is  provided  in  Figure  14. 

Because  the  laser  sheet  can  only  illuminate  the  flow  above  the  leading  edge  (a  shadow  rc'gion 
is  created  below  the  leading  edge),  the  trailing  edge  jet  is  emitted  from  the  lower  blowing  slot 
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so  that  any  stagnation  point  inovcmeiit  occurs  on  the  upper  surface  of  the  leading  edge.  The 
flow  passing  over  the  slot  lip  is  of  interest  as  well,  so  tests  are  repeated  with  upper  slot  blowing, 
since  the  light  sheet  also  illuiiiiiiates  the  flow  above  the  upper  surface  lip.  Thus,  since  the 
PIV  experiments  utilize  both  upper  and  lower  slot  blowing,  it  is  important  to  compare  Cp 
distributions  for  both. 

Figure  15  compares  Cp  distributions  for  upper  and  lower  slot  blowing  in  the  closed  test 
section  configurations.  Surprisingly,  there  are  substantial  differences  in  the  Cp  distributions 
for  upper  and  lower  slot  blowing  when  the  te.st  section  is  enclosed  for  PIV  testing.  With 
upper  slot  blowing  (blowing  on  the  same  side  as  the  j)olycarboriate  ceiling),  the  Cp  values  are 
more  negative.  Figure  15  also  reveals  that  the  Cp  data  for  lower  slot  blowing  in  the  PIV  test 
configuration  (blowing  on  the  .same  side  as  the  foam  floor)  agree  with  the  upper  slot  blowing  Cp 
data  obtained  when  the  test  section  is  fully  enclosed  with  foam.  Lower  slot  blowing  experiments 
are  not  performed  with  the  all-foam  closed  test  section. 

To  rule  out  the  blowing  slots  as  the  cause  for  these  differences,  Cp  distributions  for  upper  and 
lower  slot  blowing  in  the  open  test  section  are  compared  in  Figure  16  for  similar  test  conditions. 
The  slight  differences  between  the  two  distributions  are  minor,  and  it  can  be  concluded  that  ihe 
variation  in  the  Cp  distributions  in  Figure  15  is  not  the  result  of  dissimilar  blowing  between  the 
upper  and  lower  slots.  Instead,  it  appears  that  the  boundary  conditions  imposed  by  the  tunnel 
surface  on  the  same  sid('  as  the  blowing  slot  have  a  considerable  influence  on  the  flow  field.  A 
potential  flow  analysis  is  performed  to  provide  insight  into  this  observed  flow  behavior. 


4.1.4  Potential  Flow  Analysis 

A  full  description  of  the  potential  flow  analysis,  including  MATLAB  codes,  is  provided  in 
Wetzel  (2011).  The  approach  is  outlined  in  this  section.  Potential  lifting  flow  over  a  cylinder  in 
a  frees treain  can  be  transformed  to  flow  over  an  ellipse  via  conformaJ  mapping  (Pant on  2005: 
Katz  k,  Plotkin  2001).  Consider  the  z  =  x  iy  plane  whose  origin  is  located  at  the  center  of  a 
cylinder  of  ratlins  R.  The  7/  plane,  whose  origin  is  located  at  the  center  of  the  ellipse,  is  defined 
by 

o 

n  =  (11) 

where  a  is  the  transformation  constant  ^^Titten  in  terms  of  the  the  cylinder  radius  /?,  the  ellipse 
semi-major  axis  A,  and  the  ellipse  semi-minor  axis  B. 

n  =  \{A  +  D)  (12) 

a=^^R{A-B)  (13) 


The  complex  potential  for  lifting  flow  around  a  cylinder  of  radius  R  is  the  suininatioii  of  the 
j)otentials  for  a  uiiiforni  freestream,  a  doublet,  and  a  vortex. 


F{z)  = 


4- 


iV  111  Z 
27r 


(14) 
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Taking  the  derivative  of  Equation  14  yields  the  complex  velocity. 


(IF 


•f 


/T 

2t:z 


(15) 


To  determine  the  surface  pressure  on  an  ellipse  with  circulation  in  a  freestreain,  the  coini^lex 
velocity  W{z)  must  be  transformed  to  1T(7;)  using  the  conformal  mapping  function  given  by 
Equation  11.  Thus,  it  follows  that. 


\V{r^)  = 


drj 


(i^(h 
dz  di] 


(K)} 


and  tlie  complex  velocity  on  the  ellipse  surface,  where  ;  =  is  given  by 


ir(r;)  = 


Uoci 


,i2W 


1)  + 


pi29  _  ai 


(17) 


Since  li^(t;)  =  u  —  iv,  the  magnitude  of  the  velocity  on  the  surface  and  hence  the  C',,  distribution 
along  the  surface  of  the  ellipse  can  be  computed. 


Cp  =  1  - 


t/2  +  1,2 

172 

DC 


(18) 


The  aforementioned  analysis  can  be  extended  to  include  rigid  boundaries,  like  ground  and 
coiling  pianos,  using  the  nictliod  of  image's  (Katz  &  Plotkiii  2001).  Consider  potc’iitial  lifting  flow 
around  a  circular  cylinder  placed  a  distance  he  from  a  ceiling  plane.  An  ’image"  cylinder  with 
circulation  of  opposite-sense  from  the  real  cylinder  is  ’’placed”  a  distance  he  from  the  opposite 
side  of  the  plane.  The  origin  of  the  real  cylinder  coincides  with  the  origin  of  the  Zi  =  xi  iyi 
plane,  and  the  origin  of  the  image  cylinder  coincides  with  the  origin  of  the  22  =  ^2  +  ?y2  plane, 
which  is  offset  from  the  zi  plane  such  that  zo  —  z\  —  i2he.  The  mapping  transformation  for 
the  real  cylinder  is  /71  =  21  H~  a^fz\,  and  the  mapping  transformation  for  the  image  cylinder  is 
rj2  =  2*2  +  fz2-  The  complex  potentials  for  the  real  and  imaginary  cylinders  are  given  by  the 
following  equations. 


Fi(z,)  =  U^ 


+ - 1., 


f2(=2) 


UpcR^ 

2*2 


(19) 

(20) 


An  expression  is  desired  for  the  complex  velocity  of  the  flow  around  both  ellipses  as  a  fuiK'tion 
of  r/i.  For  the  real  ellipse,  this  is  straight for\vard  and  follows  the  steps  outlined  in  the  freestreain 
flow  analysis. 


dFi  dzi 


+ 


?r  ~ 

27r2i 


(21) 


The  complex  velocity  for  the  image  ellipse  is 


df-2  dz2  drf2 
dZ2  (17/2  (77/1 


(22) 
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The  distance  between  the  cylinders  and  the  distance  between  the  ellipses  are  not  the  same,  since 


(23) 


Therefore,  770  =  ^/i  “  2ihp,  dr/o/dz/i  =  1,  and  the  complex  velocity  for  the  image  ellipse  is  given 

by 


,  iV  \  1 

^2  27r22  /  1  -  2^ 


(2}) 


Substituting  Z2  zi  -  2ihc  and  then  writing  both  complex  velocities  in  terms  of  Zi  =  /?r*^ 
yields  the  following  complex  velocities. 


W2{rn)  = 


Vri(77i)  = 

1 

1  ~ 

^  {Re^^-2tkcy 


^oc  -  1)  + 


pi20  _ 

e  ^ 

U^R- 


+ 


iV 


{Re^f>  -  2f/io)2  2TTiRe'”  -  2iK) 


(25) 

(2(i) 


Finally,  the  complex  velocity  on  the  surface  of  the  real  ellipse  is  just  the  superposition  of  tlu* 
image  and  real  ellipse  complex  velocities, 


ir(r,i)  =  W^i{t/i)  +  U2{r/i). 


(27) 


Cp  is  again  computed  using  equation  18. 

If  a  ground  plane  is  present  instead  of  a  ceiling  plane,  then  the  only  change  to  the  previous 
analysis  is  that  Z2  =  Z\  2ihc  and  7/2  =  771  +  2ific.  The  presence  of  both  ground  and  ceiling 
planes  is  accounted  for  by  an  infinite  nuniber  of  image  ellipse's,  since  each  image  itself  is  reflc'cted 
by  the  additional  plane  (Katz  h  Plotkin  2001).  Thus,  the  complex  velocity  on  the  surface  of 
the  real  ellipse  between  two  planes  is  given  by  the  following  summation. 


C.O  (e'”  -  1)  + 
1 

00  - 


»'{rn)  = 


U^R^ 


+ 


{l)n+llT 


?;r 


n=l 


(Re^^  —  2r\ihcY  27r(/?c*^  —  2nihc) 
+  27177ic)^  2i:{Re^^  +  2mhc)  _ 


(28) 


For  moderate  circulation,  the  lift  coefficient  computed  from  the  results  of  Ecpiation  28  converges 
to  within  three  significant  digits  with  well  under  500  image  ellipses  considered.  However,  this 
result  is  largely  dependent  on  the  specified  circulation,  as  convergence  is  established  with  few 
as  20  images  for  certain  lift  values. 

Cp  distributions  for  potential  lifting  flow  around  a  20%  ellipse  in  a  variety  of  flow  configura¬ 
tions  are  compared  in  Figure  17  for  a  freestream  c/  =  1.90.  Specifically,  the  impact  of  a  ground 
plane,  ceiling  plane,  and  both  ground  and  ceiling  planes  is  assessed.  The  ground  and  ceiling 
planes  are  placed  at  the  same  distance  from  the  ellipse,  whose  chord  is  identical  to  the  airfoil 
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used  ill  experiment s,  as  the  floor  and  ceiling  installed  in  the  UFAFF.  Figure  17  reveals  a  few 
interesting  details.  First,  the  addition  of  a  ceiling  plane  drastically  augments  the  lift,  whether 
or  not  a  ground  plane  is  present.  Rirthermore,  the  addition  of  a  ground  plane  alone  actually 
reduces  the  lift  compared  to  the  freestream  case. 

The  potential  flow  analysis  only  assumes  impermeable,  rigid  boundaries,  whereas  in  actual 
testing,  a  combination  of  impermeable  (polycarbonate)  and  porous  (foam)  boundaries  are  used. 
For  this  reason,  while  the  potential  flow  results  cannot  be  expected  to  match  experiments 
exactly,  the  trends  revealed  by  Figure  17  can  still  be  applied  to  the  measured  observations.  P'or 
example,  Figure  17  indicates  that  the  presence  of  an  impermeable  ceiling  plane  produces  the 
greatest  suction-side  and  lowest  pressure-side  Cp  magnitudes,  regardless  of  whether  a  ground 
plane  is  present,  absent,  or,  hence,  impermeable  or  porous.  This  agrees  with  the  measured  data 
presented  in  Figure  15,  where  the  greatest  suction-side  and  lowest  pressure-side  Cp  magnitudes 
are  observed  when  blowing  is  utilized  on  the  same  side  as  the  polycarbonate  tunnel  ceiling. 
Additionally,  Figure  15  shows  little  variation  in  the  measured  Cp  distributions  when  blowing  is 
utilized  on  the  same  side  as  the  porous  foam  ceiling  and  the  pressure-side  boundary  is  either 
polycarbonate  or  foam.  This  agrees  with  the  theory  in  Figure  17,  which  indicates  that  the 
presence  of  a  ground  plane  has  a  minimal  influence  on  the  Cp  distribution. 

To  further  illustrate  these  tix'iids,  the  potential  flow  solutions  for  the  ellipse  near  a  ground 
plane  and  ceiling  plane  alone  are  compared  with  measiirenients.  Since  the  potential  flow  .solu¬ 
tions  do  not  account  for  the  pres<mce  of  the  trailing  edge  jet  and  its  contribut  ion  to  the  trailing 
edge  suction  peak,  measured  and  computed  Cp  distributions  cannot  be  compared  using  equiv¬ 
alent  lift  coefficients.  Instead,  the  theoretical  and  measured  distributions  are  match<'d  at  the 
leading  edge  stagnation  points.  The  measured  closed  test  section  Cp  distributions  for  blowing  on 
the  same  side  as  and  opposite  side  from  the  polycarbonate  tunnel  ceiling  are  compared  to  tlu'ory 
for  flow  with  a  ceiling  plane  and  ground  plane,  respectively,  in  Figure  18.  The  measured  leading 
edge  stagnation  points  axe  found  to  match  with  theory  for  both  data  sets  when  the  fn^estreain 
lift  coefficient  is  1.90.  At  the  leading  edge,  the  theory  agrees  well  with  the  measurements,  and 
the  trends  revealed  by  Figure  17  are  clearly  followed  by  the  measured  data. 

Further  potential  flow  analysis  indicates  the  difference's  between  the  computed  closed  te'st 
section  and  freestream  surface  pressure  distributions  are  insignificant  when  the  tunnel  blocloige 
ratio,  computed  from  the  ratio  of  the  airfoil  thickness  to  tunnel  height,  is  less  than  6%.  By  com¬ 
parison,  the  tunnel  blockage  ratio  for  the  present  investigation  is  14%,  and  the  blockage  ratios 
for  studies  by  Abramson  (1975)  and  Novak  &  Cornelius  (1986)  were  8%  and  7%,  respectiv('ly. 

4.1.5  PIV  Results 

Although  theory  and  measurements  agree  that  there  are  large  overall  differcmces  in  the  Cp 
distributions  depending  on  the  suction-side  boundary  condition,  in  general,  enclosing  the  test 
section  clearly  produces  the  leading  edge  suction  peak  absent  in  the  open  test  section  mea¬ 
surements.  The  influence  of  ('iiclosing  the  test  section  on  the  leading  edge  flow  is  captured  by 
PIV.  Mean  flow-  streamlines  for  both  open  and  closed  test  section  configurations  at  a  variety 
of  lower  blowing  slot  momentum  coefficients  are  shown  in  Figure  19.  In  the  open  test  section 
configuration,  leading  edge  stagnation  point  movement  and  streamline  curvature  are  minimal 
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as  increases.  Conversely,  when  the  test  section  is  enclosed,  the  leading  edge  stagnation  point 
moves  considerably  with  increasing  and  streamline  curvature  is  greatly  enhanced. 

The  primary  motivation  for  this  research  investigation  is  acoustics.  It  has  been  revealed 
that,  upstream  of  the  slot,  the  Cp  distributions  measured  with  the  test  section  opt'ned  and 
clo.sed  differ  significantly.  Passive-slot  noise  is  theorized  by  Howe  (2002)  to  be  produced  by 
turbulence  in  the  freest  ream  boundary  layer  .scattering  off  the  slot  lip.  If  the  characteristics  of 
the  turbulent  boundary  layer  near  the  slot  vary  with  test  section  configuration,  then  so  too  does 
the  sound  produced.  The  turbulent  boundary  layer  passing  over  the  slot  lip  is  measured  using 
PIV  in  both  open  and  closed  tunnel  configurations  for  a  variety  of  upper  slot  blowing  momentum 
coefficients.  The  results  are  presented  in  Figures  20,  21,  and  22.  Profiles  of  mean  tangential 
and  normal  velocities,  turbulence  intensities,  and  Reynolds  stress  are  plotK'd.  For  clarity,  only 
every  fourth  datum  point  is  displayed,  and  uncertainty  bounds  are  represented  by  solid  lines. 
Without  blowing  (Figure  20),  the  profiles  are  normalized  using  the  local  freest  ream  velocity  and 
the  boundary  layer  thickness.  With  blowing  (Figures  21  and  22),  the  mean  tangential  velocity 
profiles  appear  more  like  wall  jet  profiles,  and  the  local  maximum  velocity  and  the  boundary 
layer  thickness  based  on  where  the  velocity  is  0.99f/„,ax  QJ'e  used  to  normalize  the  data.  For 
all  cases,  there  are  iiiinor  differences  between  the  profiles,  and  these  deviations  are  typically 
within  the  iiieasurenieiit  uiicertaiiity.  From  an  acoustics  perspective,  tiie  differences  between 
the  profiles  should  be  negligible,  and  the  sound  produced  by  the  interaction  of  the  turbulent 
boundary  flow  with  the  slot  lip  should  be  the  same  whether  or  not  the  test  section  is  opened  or 
closed. 

4.2  Curved  Wall  Jet  Flow 

The  primary  mechanism  exploited  by  circulation  control  is  the  curved  wall  jet,  which  is  ('xpected 
to  play  an  important  role  in  the  production  of  circulation  control  noise.  Curvature  noise  is 
theorized  by  Howe  (2002)  to  be  produced  by  the  interaction  of  boundary  layer  turbulence  with 
the  rounded  trailing  edge,  and  the  scales  required  to  assess  Howe’s  model  of  this  noise  -  the 
displacement  thickness,  mean  velocity,  and  friction  velocity  -  are  clearly  going  to  be  dictated 
by  the  interaction  of  the  curved  wall  jet  with  the  freestream.  Thus,  there  is  a  need  for  liigh- 
resolution  flow'  measurements  of  the  curved  w'all  jet  to  determine  these  s(*ales.  Of  greater 
importance  and  interest,  however,  is  the  potential  for  flow  similarity.  If  the  circulation  control 
wall  jet  exhibits  similarity,  then  it  may  be  possible  to  predict  the  evolution  of  the  flow'  and 
its  various  length  and  velocity  scales.  Such  a  finding  w'ould  have  direct  application  to  Howe's 
model  and  also  the  design  of  circulation  control  airfoils  for  a  variety  of  applications. 

Unfortunately,  the  literature  provides  evidence  that  full  flow  similarity  may  only  be  possible 
for  very  specific  geometries.  Guitton  Newman  (1977)  presented  a  similarity  solution  fur  a 
curved  w'all  jet  in  the  absence  of  a  freestream  and  concluded  that  full  flow  similarity  is  only 
possible  if  the  surface  is  defined  by  a  logaritliniic  spiral,  c.g.  r  cx  where  R  is  the  local 

radius  of  ciirvature,  and  x  is  the  arc  length  from  the  origin.  A  logarithmic  spiral  is  illustrated 
in  Figure  23.  Their  analysis  is  extended  in  the  next  section  for  the  curved  w^all  jet  with  an 
external  freestream. 
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4*2.1  General  Similarity  Solution 


Using  the  turbulent  forms  of  the  governing  equations  in  polar  coordinates,  Guitton  &  Newman 
(1977)  developed  a  similarity  solution  for  a  curved  wall  jet  in  the  absence  of  a  freestn^am, 
assuming  the  following  forms  of  the  similarity  functions  and  512- 


W  —  UmaxS  iv) 

(2!)) 

-U'v'  =  U^nax9n{l) 

(30) 

u'2  _ 

(31) 

y 

7] 

(32) 

Umax  the  local  maximum  velocity,  and  y\/2m  Is  the  local  normal  distance  from  tin*  surface 
where  the  velocity  is  ^Umax-  Guitton  and  Newman  showed  that  the  wall  jet  flow  can  only  be 
self-similar  provided  the  surface  is  defined  by  a  logarithmic  spiral,  r  ex  The  ratio  .1  /R 

dictates  the  rate  at  which  the  local  radius  of  curvature,  R,  grows.  The  larger  the  value  of  x/R, 
the  larger  the  local  radius  of  curvature  at  some  arc  length  x.  Table  7  lists  the  relevant  length 
and  velocity  scales  and  how  each  scales  with  arc  length  x. 

With  an  external  freestream,  the  velocity  may  never  decay  to  ^Umax-  hor  such  conditions. 
Launder  &  Rodi  (1983)  suggested  the  use  of  a  defect  velocity,  defined  as  Umax  —  Ue,  where  L\ 
is  the  velocity  where  the  Reynolds  stress  decays  to  a  negligible  value,  and  the  length  y^.1/2, 
defined  as  the  normal  distance  from  the  .surface  where  the  velocity  is  ^{Umax  d-  Ue)^  These 
scales  are  illustrated  in  Figure  24.  The  analysis  of  Guitton  &  Newman  (1977)  is  extended  for  a 
curved  wall  jet  in  the  presence  of  a  freestream  by  assuming  the  following  arbitrary  forms  of  the 
similarity  functions. 


»  -  f/,/'(r;)  +  t/o 

(33) 

-xi'v'  =  Ufgi2{v) 
j/'2  _  1)12  =  Ulg{g) 

(34) 

(3.5) 

y  -  92 

1]  = 

y\ 

(36) 

If  the  scales  suggested  by  Launder  k,  R.odi  (1983)  are  applied  to  similarity  functions  listed  in 
Equations  33  to  36,  then  U\  =  Umax  “  Up,  Uo  =  Up,  y\  —  —  ymax^  ^nd  yy  =  ymax- 

The  remainder  of  the  similarity  analysis  follows  the  same  steps  outlined  by  Guitton 
Newman  (1977)  and  is  included  in  detail  in  Wetzel  (2011).  The  similcirity  equation  is 


fi2 

t/i  iix 


wi  dUo  .  ...  dui  U')  .  ...  dt/o  [Uo 


fh2 


yi  dUo  f  .ff.  dy\  L2  f  ,  ^y2 

U,  t/,  dx 


~  ■^(512)  —  0 


(37) 
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where 


4 .  (j'j . .  /  V«)  4  (/, .  /  V..) , .  /■  (?|)\  1.« 


B 


C  =  ivf)  + 

^  =  ^if)  - 

Ui^-’ ’  ul 


(40) 

(41) 


When  i/2  =  y2  =  0,  Equation  37  reduces  to  the  form  of  the  similarity  solution  provided  by 
Guitton  &  Newsman  (1977).  In  the  general  case,  all  terms  outside  the  parentheses  in  Equation 
37  must  be  constant.  Consider  the  term, 


^(r/'+W). 


;-i2) 


Since  dyxjdx  must  be  constant,  it  follows  that  y\  oc  x.  Similarly,  note  the  term, 


^(wi2  +  2</i2).  (43) 

Since  y\/R  must  be  constant,  then  R  (x  y\  (x  x.  Therefore,  in  order  for  true  .similarity  to  be 
achieved,  the  surface  must  be  a  logarithmic  spiral,  i.e.  r  x  where  R  is  the  local  radius  of 

curvature  labeled  in  Figure  23.  Table  7  lists  the  length  and  velocity  scales  from  this  analysis 
and  how  each  scales  with  arc  length  x. 

Since  the  trailing  edge  of  the  circulation  control  airfoil  has  a  constant  radius,  complete  self¬ 
similarity  of  the  flow  cannot  be  achieved  in  the  present  investigation.  A  log-spiral  trailing  edge 
is  not  practical  for  underwater  vehicle  applications  where  control  surface  symmetry  is  desired. 
In  addition,  a  log-spiral  design  would  likely  interfere  with  the  secondary  blowing  slot  used  to 
control  excessive  jet  attachment  at  extremely  high  primary  blowing  slot  inoinentum  coefficients. 

Although  full  flow  similarity  is  not  possible  in  the  present  investigation,  even  partial  flow 
siinilarity  could  still  prove  useful.  For  instance,  similarity  of  the  mean  tangential  velocity  alone 
would  potentially  provide  a  means  of  estimating  the  displacement  thickness,  mean  velocity, 
and  perhaps  friction  velocity  required  for  Howe’s  model.  Before  the  curved  wall  jet  flow  is 
examined  for  similarity,  dimensional  analysis  is  used  to  identify  the  dependent  parameters  in 
this  investigation. 


4.2.2  Dimensional  Analysis 

A  dimensional  analysis  is  performed  to  determine  the  important  dimensionless  parameters  as¬ 
sociated  with  the  trailing  edge  flow  field.  As  an  example,  the  location  of  maxinuiin  velocity 
ymax  tised  as  the  relevant  dependent  length  scale.  Any  other  dependent  length  or  velocity 
scale  could  be  used  in  its  place. 

First,  all  significant  parameters  in  the  problem  are  identified.  Note  that  the  flow  is  assumed 
to  be  incompressible  since  Mjet  <  0.3  for  all  cases  considered.  The  geometric  properties  include 
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the  chord  c,  slot  height  /i,  trailing  edge  radius  and  lip  thickness  1.  Of  these,  the  chord, 
trailing  edge  radius,  and  lip  thickness  are  fixed  in  the  curr<*nt  experinu'iital  investigation.  Flow 
properties  include  density  p,  viscosity  /i,  freestream  velocity  Uooy  velocity  Ujet,  and  //moj-* 

If  c,  Uoo^  and  p  are  chosen  as  the  repeating  parameters,  the  II  groups  as.semblod  are:  111  = 
/i/c,  rio  =  r/c,  ria  =  l/c.  U4  =  pU^clp  =  Rccs  Us  =  Ujet/U^,  and  Uq  =  ymax/<^’  II5  can  be 
recast  into  the  momentum  coefficient, 


n;  =  2n,(n,)^  =  2P)l^ 


(.H) 


Note  that  the  jet  Reynolds  number,  defined  as 

Re,,,  =  (.15) 

P 

IS  related  to  the  chord  Reynolds  number,  inonientuni  coefficient,  and  slot  heiglit-to-cliord  ratio. 


Hence,  of  and  /i/c,  only  three  are  independent.  Therefore,  the  dimensionless 

length  scales,  such  as  ymax/c  and  velocity  scales,  like  Umax/^x^  arc  dependent  on  the  following 
parameters. 


c  \c  c  c 


(17) 


As  noted  earlier,  r/c  and  l/c  are  constant  in  the  present  investigation.  However,  /i/c, 
and  Rcc  are  variable.  Thus,  two  of  these  three  parameters  can  be  kept  constant,  and  the  effect 
of  varying  the  third  parameter  on  ymax  Q^nd  other  length  and  velocity  scales  can  be  studied. 

Table  3  summarizes  the  six  different  cases  that  are  included  in  this  study  of  the  curved  wall 
jet.  Cases  one  through  three  represent  tests  where  Rcc  and  h/c  are  fixed,  but  changes.  In 

cases  four  and  five,  and  h/c  are  fixed,  but  Rcc  varies.  Finally,  in  cases  four  and  six, 

and  Rcc  are  fixed,  but  h/c  differs.  In  the  following  sections,  results  from  these  six  cases  are 
presented  and  discussed. 


4.2.3  Flow  Characteristics 

Before  studying  the  PIV  measurements  of  the  curved  wall  jet,  the  surface  pressure  along  the 
Coanda  surface  is  examined.  Trailing  edge  Cp  data  for  the  cases  listed  in  Table  3  are  plotted 
in  Figure  25.  Initially,  a  favorable  pressure  gradient  accelerates  the  flow  for  all  cases  except 
=  0.0039,  where  Ujet  ^  Shortly  thereafter,  near  rO/c  =  0.02,  a  .strong  adverse  ]>ressure 
gradient  develops.  The  trends  displayed  in  Figure  25  follow  similar  findings  reported  by  Nov’ak 
et  al.  (1987).  Trailing  edge  Cp  provides  insight  into  the  expected  similarity,  or  lack  thereof, 
of  the  curved  wall  jet  flow.  In  their  study  of  a  curved  wall  jet  in  the  absence  of  a  freestream, 
Neuendorf  Wygnaiiski  (1999)  observed  full  similarity  of  0  (inner  and  outer  regions)  in  the 
constant  pressure  region  but  only  outer  region  similarity  of  U  in  the  adverse  pre.ssnre  gradient 
region. 
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Profiles  of  iiioaii  velocity  {0  and  V),  turbulence  intensity  and  and  Reynolds 

stress  from  PIV  measurements  along  the  rounded  trailing  edge  are  plottcxl  in  Figure  26 

at  a  variety  of  downstream  distances  from  the  slot  for  Her  =  6.5 x  10^,  h/c  =  0.0019,  and  — 
0.0039,  0.015,  and  0.057.  Profiles  are  not  included  for  =  0.0039  after  rO/c  =  0.040  because 
separation  occurs  upstream  of  that  position.  The  distance  from  the  surhice  is  iioruialized  using 
the  slot  height  h.  and  the  flow  quantities  are  normalized  using  the  jet  velocity  Uj^t-  Uncertainty 
bounds  are  removed  from  the  plots  for  clarity.  The  mean  tangential  velocity  profiles  reveal  how 
the  jet  decays  and  spreads  as  it  travels  away  from  the  slot.  As  U  decays,  the  mean  normal  ve¬ 
locity  V  increases  and  separation  is  approached.  The  turbulence  intensities  and  Reynolds  stress 
are  significant  in  the  iniinediate  vicinity  of  the  blunt  slot  lip  and  the  high-shear  region  that 
exists  there.  As  the  jet  spreads  and  entrains  fluid  from  the  freestream,  Reynolds  stress  becomes 
non-negligible  over  a  growing  region  of  the  trailing  edge.  The  similarity  of  the  measured  curved 
wall  jet  flow  is  assessed  in  the  following  section. 

4.2.4  Flow  Similarity 

Novaik  Cornelius  (1986)  and  Novak  et  al.  (1987)  showed  that  the  outer  region  of  U  along  the 
trailing  edge  of  their  circulation  control  airfoil  exhibited  similarity  when  normalized  using  the 
scales  suggested  by  Launder  &  Rodi  (1983)  (refer  to  Figure  24).  These  same  scales  are  applied 
to  the  current  data  set  in  order  to  determine  if  similarity  is  achieved.  In  addition,  the  data  is 
normalized  using  the  scales  found  by  Zhou  &  Wygnaiiski  (1993)  to  collapse  the  outer  region  of 
U  for  a  plane  wall  jet  in  an  external  freestream.  Those  scales  include  local  Umas^  t'oot  Vinaxy 
and  which  is  the  normal  distance  where  U  =  In  the  current  set  of  experiments, 

Uoc  is  replaced  with  f/mm,  the  velocity  where  the  outer  region  mean  shear  approaches  zero,  and 
i^  replaced  with  y,n,i/2^  the  normal  distance  where  U  =  ^{Umax  since  the  velocity 

does  not  always  decay  to  values  of  ^Umax-  These  scales  are  illastrated  in  Figure  27. 

To  normalize  the  flow  profiles,  the  various  length  and  velocity  scales  described  must  first  be 
determined.  Umax  is  found  by  detecting  the  point  of  maximum  velocity  in  each  profile,  fitting 
a  Gaussian  curve  to  this  point  and  its  adjacent  four  neighbors,  and  then  finding  the  maximnm 
of  the  Gaussian  curve.  With  Umax  known,  //mar  is  dcteriniiied.  To  find  ymin^  mean  shear 
{ilU /<ly)  is  computed  starting  at  ymax  and  moving  outwards  from  the  surface  until  the  sign  of 
the  shear  flips,  and  then  a  third-order  polynomial  is  fit  to  the  two  points  on  either  side  of  the 
zero-crossing,  and  the  exact  zero-crossing  is  determined  and  set  as  ymin-  For  a  f('w  profiles, 
the  shear  never  changes  sign  in  the  outer  region.  In  those  situations,  ymin  i^  specified  as  the 
first  instance  where  {h/Ujf,t){dU /dy)  <  4  x  10”®  (or  U fUj^t  changes  by  less  than  6  x  10“^). 
This  particular  shear  threshold  is  determined  by  inspection  of  the  data.  With  ymm  sptrified, 
Umtn  is  found  by  evaluating  a  third-order  polynomial  fit  of  the  velocity  data  at  ymm-  Fm.i/2 
computed  once  Umax  ^nd  Umin  are  known,  and  then  t/m,i/2  determined  by  a  piecewise  cubic 
interpolation  along  all  the  points  between  ymax  and  ymin^  To  determine  t/r,  first,  is 

determined  in  similar  fashion  to  Umax  using  a  Gaussian  curve  fit.  Then  the  first  instance  where 
n'v* /(u'v')max  <  0.05  (this  value  is  also  chosen  by  inspection;  slightly  higher  thresholds  do  not 
significantly  affect  the  results)  above  the  position  of  maximum  Reynolds  stress  is  deemed  the 
negligible  Reynolds  stress  location.  A  third-order  polynomial  fit  of  the  velocity  data  is  then 
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evaluated  at  this  location  to  determine  Ue.  ^eA/2  follows  the  same  approach  used  to  find  (/n».i/2* 

The  length  scales  ymax^  J/mim  l/e^  ymA/2^  j/e,i/2  are  plotted  with  the  mean  tangential 
velocity  and  Reynolds  stress  profiles  for  =  0.015  in  Figure  28.  Profiles  are  provided  from 
rO/c  =  0.0027  until  just  past  separation  at  rO/c  =  0.047.  As  Figure  28  attests,  there  is  little 
differonre  between  ymin  and  and  ymA/2  and  yej/2- 

The  scales  are  used  to  normalize  the  outer  regions  of  the  mean  tangential  velocity  profiles  in 
Figure  29,  where  it  is  seen  that,  using  either  the  scaling  suggested  by  Launder  &  Rodi  (1983)  or 
the  modification  to  the  approach  by  Zhou  <L'  Wygnanski  (1993),  the  outer  streamwise  profiles 
exhibit  similarity.  Note  that  the  velocity  data  at  rO/c  =  0.0027  are  not  included,  as  they  do 
not  collapse  with  the  other  profiler.  In  fact,  only  profiles  in  the  constant  pressure  or  adverse 
pressure  gradient  region  along  the  trailing  edge  exhibit  similarity,  agreeing  with  the  observations 
of  Neuendorf  &;  Wygnanski  (1999),  Novak  &  Cornelius  (1986),  and  Novak  et  al.  (1987).  The 
solid  line  in  Figure  29  represents  a  hyperbolic  tangent  of  the  form 

[7=0.5  [l  -  1.315 tanh(?/'l)]  ,  (48) 

where  U  and  y  <ire  the  normalized  velocity  and  distance  from  the  surface,  respectively. 

Now  that  the  outer  region  of  U  has  been  found  to  exhibit  similarity,  the  inner  region  of  U 
is  studied.  It  is  difficult  to  measure  extremely  close  to  the  surface  with  PIV  (recall,  the  slot 
height  is  typically  1  inin).  so  inner  region  data  are  limited  compared  to  the  outer  region  of  the 
flow.  Nevertheless,  there  are  sufficient  data  to  deteniiiue  whether  or  not  siiuilarity  is  iuhieved. 
Inner  region  U  data  are  normalized  using  ymax  Umax  and  plotted  in  Figure  30.  Unlike  the 
outer  region,  the  inner  region  does  not  exhibit  similarity.  These  results  agree  with  Neuendorf 
&  Wygnanski  (1999),  who  found  that  the  inner  region  of  U  for  a  curved  wall  jet  in  quiescent 
surroundings  does  not  exhibit  similarity  in  the  adverse  pressure  gradient  region. 

The  remaining  flow  profiles,  V ,  (u'^)^/^,  (e'2)^/~,  and  uV,  also  do  not  exhibit  similarity 
using  any  eombination  of  mean  velocity  and  length  scales.  However,  the  similarity  of  U  in  the 
outer  region  alone  may  permit  the  prediction  of  that  portion  of  the  flow  and  subsequently  the 
mean  tangential  velocity  scales,  which  are  required  to  assess  Howe’s  model  (Howe  2002).  In 
addition,  knowledge  of  ymax  1^^  useful  in  providing  at  least  a  rudimentary  predictit)ii  of  the 
local  displacement  thickness.  In  the  following  section,  the  length  and  velocity  scales  found  to 
collapse  the  outer  region  of  U  are  analyzed  for  the  six  cases  listed  in  Table  3. 

4.2.5  Length  and  Velocity  Scales 

The  set  of  length  and  velocity  scales  ba.sod  on  the  mean  shear  found  to  sufficiently  collap.se  tlu' 
outer  streamwise  profiles  do  not  require  knowledge  of  the  flow  turbulence  and  an'  thus  more 
practical  in  application.  As  such,  they  are  the  focus  of  the  remainder  of  the  analysis.  The  length 
and  velocity  scales  associated  with  this  approach,  ymax^  ?/m,i/2’  Umax^  Umin^  plotted 
ill  increments  of  r9/c  =  0.002  for  —  0.015  in  Figure  31.  The  data  are  plotted  from  just 
downstream  of  the  slot  exit  to  separation,  detected  by  the  first  sign  of  flow  reversal.  Figun'  31 
further  illustrates  the  development  and  spread  of  the  jet.  The  velocity  scales  initially  increase, 
reach  a  maximum  value  near  r9/c  =  0.015,  then  decay  at  similar  rates.  The  length  scales  ymax 
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and  t/rn, 1/2  nearly  constant  initially,  then  iiicreiise  gradually,  and  finally  grow  sharply  as 

separation  is  approached.  On  the  other  hand,  ymin  increases  immediately  and  then  follows  a 
parabolic  form,  similar  to  the  other  length  scales. 

For  all  test  cases  listed  in  Table  3,  the  trends  of  t/rnax?  i/m, 1/2^  Umax^  Umin  anal>^ed  in 
an  attempt  to  collapse  each  set  of  curves.  Only  data  where  the  outer  region  exhibits  similarity 
is  included.  Neuendorf  &  Wygnanski  (1999)  found  that,  for  a  curved  wall  jet  in  the  absence 
of  a  freestream,  the  decay  of  Umax  ai^d  the  rate  of  spread  of  the  jet  described  by  the  location 
where  U  =  \Umax  could  be  scaled  using  the  local  kinematic  jet  momentum, 

j  =  2/m/2  r  ( 7j^)  d  (^)  =  0.78t/^„^  (19) 

and  the  wall  radius  in  the  constant  pressure  region.  However,  unlike  the  curved  wall  jet  of 
Neuendorf  and  Wygnanski,  the  entire  velocity  profile  of  the  circulation  control  curved  w'all  jet 
does  not  exhibit  similarity.  Hence,  the  local  kinematic  jet  momentum  cannot  be  exptx'ted  to 
collapse  the  scales.  Likewise,  Neuendorf  and  Wygnanski  observed  that  the  rate  of  spread  could 
be  expressed  solely  as  a  function  of  the  wall  radius  and  downstream  angular  position  from  t  he 
slot.  This  is  clearly  not  observed  for  the  circulation  control  airfoil,  as  shown  in  Figures  32  and 
33.  For  a  planar  w'all  jot  in  an  external  freestream,  Zhou  k.  Wygnanski  (1993)  found  that  the 
length  and  velocity  scales  describing  the  planar  wall  jet's  spread  and  decay  could  be  expressed 
as  a  function  of  a  normalized  downstream  distance  from  the  slot  defined  by  where  A'  is 

the  dimensional  distance  from  the  slot,  u  is  the  kinematic  viscosity,  and  J  =  h{Ujet  —  UocWjtt 
is  the  excess  of  kinematic  momentum  flux  near  the  nozzle.  A  dimensionless  velocity  ratio 
{Ujet  —  Ucyo)/{Ujf>f  Uoo)  was  also  used.  These  parameters  do  not  collapse  the  current  data  set 
either.  Instead,  the  data  are  best  fit  using  Rec,  Cfi,  and  h/c.  as  shown  in  Figures  32  through 
35.  A  power  law  curve  is  fit  to  each  data  set  by  niiniinizing  the  square  of  the  error.  The 
scales,  as  a  function  of  rO/c,  are  found  to  collapse  with  the  product  of  Rcc  and  referred 
to  as  the  Reynolds  corrected  momentum  coefficient,  which  has  recently  been  found  to  st'ale  lift 
increnieiits  for  active  control  of  airfoil  flow'  separation  (Stahiov  k  Seifert  201U).  This  scaling 
parameter  takes  into  consideration  the  freestream  boundary  layer  and  its  development  as  a 
function  of  Rxe,  Recall  Figures  21  and  22,  where  at  higher  values  of  C^,  the  mean  tangential 
velocity  profile  of  this  npst regain  '^boundary  layer”  does  not  actually  resemble  a  conventional 
boundary  layer.  The  Reynolds  corrected  momentum  coefficient  provides  .some'  basis  to  account 
for  this  effect.  Best-fit  equations  for  the  length  and  velocity  scales  are  provided  in  Equations 
50  to  53. 
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Equations  50  and  51  reveal  approximately  a  6tli  power  dependence  on  the  arc  length  r&/c. 
Contrast  that  to  the  results  of  the  general  siiihlarity  solution,  which  indicates  that,  for  the 
self-similar  flow  over  a  log  spiral  surface,  these  scales  should  hav'e  a  linear  relationship  with  the 
arc  length. 

Tlu*  product  of  the  chord  Reynolds  niimhor  and  moinentuni  cocfticieiit  can  also  bt'  rearranged 
as  the  product  of  the  jet  Reynolds  number  and  jet-to-freestreain  velocity  ratio. 


RecCf,  = 


_  ci^jeth  Ujai  _  r,r>„  ^jft 


(51) 


Thus,  in  the  case  of  the  circulation  control  airfoil,  it  appears  that  the  rate  of  dc^eay  and  spread 
of  the  jet  are  dependent  on  the  jet  Reynolds  number.  This  is  in  stark  contrast  to  the  findings 
of  Zhoii  &  Wygnanski  (1993).  who  determined  that,  as  long  as  Umax/U^c  >  2,  the  length  and 
velocity  scales  describing  the  planar  wall  jet  in  an  external  freestream  are  independent  of  th<' 
jet  Reynolds  number.  Note  that  U-max/^oc  >  2  for  some  profiles  in  cases  2  through  5  listed  in 
Table  3. 

Equations  50  to  53  provide  a  predictive  capability  for  the  length  and  velocity  scakrs  of  the 
mean  outer  region  flow.  Predicted  and  measured  scales  are  compared  in  the  next  section. 


4.2.6  Flow  Prediction 

A  comparison  of  the  predicted  and  measured  length  and  velocity  scales  is  presented  in  l"igiirc‘ 
36  for  Rcc  =  6.5  x  10^^,  =  0.015,  and  h/c  =  0.0019.  Overall,  the  agreement  between  tlu* 

measured  and  predicted  scales  is  good.  Umin  i^>  initially  over-predicted  in  the  favorable  pressure 
gradient  region,  but  recall,  only  data  in  the  adverse  pressure  gradient  region,  where  the  outer 
region  of  U  is  similar,  is  used  to  determine  Equations  50  to  53.  Despite  that  fact,  the  other 
scales  are  predicted  reasonably  well  in  the  favorable  pressure  gradient  region. 

While  the  agreement  betw'een  the  predicted  and  measured  flow  fields  in  Figure  36  is  promis¬ 
ing,  it  is  a  rather  expected  result  since,  after  all.  the  data  for  those  test  conditions  are  used  in 
determining  the  prediction  equations.  To  truly  test  the  prediction's  capabilities,  it  should  be 
compared  to  measurements  for  a  case  not  included  in  the  prior  analysis.  Thus,  the  predicted  and 
mciisured  scales  are  compared  in  Figure  37  for  Rec  =  1.3  x  10^,  =  0.014,  and  h/c  —  0.0029. 

Although  there  are  larger  differences  between  the  predicted  and  ineasured  scales  in  comparison 
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with  Figure  36,  the  overall  agreement  is  fair. 

To  this  point,  the  curved  wall  jet  flow  has  only  been  considered  prior  to  separation.  In  the 
iK’Xt  section,  flow  separation  is  discussed. 

4.2.7  Separation  and  Stability 

It  is  well  known  that  the  circulation  control  curved  wall  jet  is  responsible  for  entraining  freestreain 
fluid,  delaying  separation  and  increasing  circulation.  The  cause  of  separation  is  porliaps  less 
clear,  and  could  be  the  result  of  a  severe  adverse  pressure  gradient  (see  Figure  25),  meandering 
streaniwise  vortices,  observed  in  curved  wall  jets  absent  of  a  freestreain  by  Likhachev  et  al. 
(2001),  Neuendorf  et  al.  (2004),  and  Han  et  al.  (2006),  or  a  combination  of  influences. 

Separation  is  first  studied  in  the  present  investigation  by  using  PI\'  to  measure  the  trailing 
edge  flow  field  and  part  of  the  wake  region.  Mean  speed  contours  and  mean  velocity  vectors 
are  presented  for  =  0,  0.014,  and  0.057  in  Figure  38.  When  =  0,  a  large  wake  is  formed 
aft  of  the  trailing  edge,  and  the  frc'cstream  flow  appears  syiiniit'tric  about  the  cbordliiie.  \\*itli 
blowing,  the  flow  along  the  upper  surface  of  the  airfoil  remains  attached  and  separates  much 
farther  downstream,  as  shown  in  Figure  38B.  The  vectors  reveal  the  extent  of  flow  turning  due 
to  entrainment  from  the  jet.  As  further  increases  to  0.057  in  Figure  38C,  the  freestreain  is 
deflected  oven  more,  and  the  separation  point  is  shifted  farther  aw’ay  from  the  blowing  slot. 

Figure  38  only  provides  a  qualitative  description  of  flow  separation.  For  a  more  quantitative 
analysis,  the  curved  wall  jet  PIV  data  are  used  to  estimate  the  separation  location  by  detecting 
the  onset  of  flow  reversal  near  the  surface.  Data  along  the  trailing  edge  surfac*o  are  analyzed  in 
increments  of  0.13  mm,  which  corresponds  to  the  vector  resolution  of  the  data  set.  Also,  since  it 
is  not  possible  to  resolve  the  flow  field  very  near  the  surface  with  PIV,  the  separation  locations 
presented  can  only  be  regarded  as  estimates,  and  the  true  separation  location  likely  occurs  just 
upstream  of  the  estimated  location.  The  separation  locations  are  plotted  as  a  function  of  C^^Rec 
in  Figure  39,  which  indicates  that  the  separation  location  moves  farther  downstream  from  the 
slot  as  the  product  of  the  chord  Reynolds  number  and  momentum  coefficient  increases.  In 
particular,  as  C^Rec  initially  increases,  the  separation  distance  increases  at  a  high  rate,  but  at 
larger  values  of  CuRcc.  ^  significant  increase  in  the  product  of  the  two  parameters  yields  only  a 
small  delay  in  separation.  A  curve,  who.se  equation  is  given  by  E(|uatioii  55,  is  fit  to  this  data  so 
the  separation  distamv  can  be  estimated  for  a  given  Reynolds  corrected  momentum  coc'fficient. 
Equation  55  accurately  predicts  the  sepfiration  locations  for  both  cases  considered  in  Section 
4.2.6  with  an  error  of  less  than  2%  . 

=  0.0085  (55) 

It  should  be  noted  that  Equation  55  appears  to  be  only  applicable  to  the  airfoil  gt'onu'try 
of  the  present  investigation.  Attempts  to  predict  the  separation  locations  reported  by  Novak 
&:  Cornelius  (1986)  and  Novak  et  al.  (1987)  result  in  errors  of  just  under  50%  (or  6  =  65°). 
How'ever,  the  trailing  edge  geometry  of  the  circulation  control  airfoil  studied  by  No\^k  and  his 
colleagues  is  significantly  different  from  the  trailing  edge  of  the  circulation  control  airfoil  in 
the  present  investigation.  Furthermore,  Equation  55  is  only  accurate  for  Rcc  >  0.  \Vith  no 
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freest  real  11,  the  curved  wall  jet  flow  separates  well  downs  t  ream  of  the  blowing  slot. 

The  mechanism  causing  separation  has  garnered  much  attention  over  the  last  decade.  Recent 
work  by  Likhachev  et  ai  (2001),  Neuendorf  et  al.  (2004),  and  Han  et  al.  (2006)  has  revealed 
that  meandering  streamwise  vortices  are  likely  the  culprits  of  curved  wall  jet  flow  separation 
Floryan  (1986)  determined  the  inviscid  stability  criterion  for  boundary  layer  and  wall  jet  flows 
over  convex  surfaces,  given  by 


dU^ 

dy 


>  0. 


(56) 


Since  the  vv'all  jet  profile  is  non-monotonic,  whether  or  not  a  freest  ream  is  present,  it  automat¬ 
ically  violates  the  inviscid  stability  criterion.  In  particular,  dU^  jdy  <  0  in  the  outer  region 
of  the  jirofile  where  the  velocity  decays.  Therefore,  the  outer  region  of  the  wall  jet  flow  over 
a  convex  surface  is  potentially  unstable  and  susceptible  to  the  formation  of  "Gortler  vortices'* 
(Floryan  1986;  Saric  1994),  Floryan  also  performed  a  viscous  stability  analysis  and  provide^d  a 
neutral  stability  curve  for  the  Gbrtler  luirnber,  which  he  found  to  be  a  function  of  the  mfixinnim 
velocity,  boundary  layer  thickness,  kinematic  viscosity,  and  curved  wall  radius. 


G  = 


U 


(57) 


Floryan  (1986)  defined  the  boundary  layer  thickness  by 


6  = 


ure 

^mai 


1/2 


(58) 


The  Gortler  number  is  evaluated  as  a  function  of  downstream  distance  from  the  slot  for 
Rcc  =  6.5  X  10''^,  =  0.015,  and  h/c  ~  0.0019,  The  results  are  plotted  against  the  neu¬ 

tral  stability  curve  computed  by  Floryan  in  Figure  40  and  indicate  that  beyond  rO jc  —  0.0047. 
the  flow  is  unstable  to  disturbances  of  critical  wavelengths  indicated  by  the  region  of  the  figure 
labeled  ‘'unstable,'*  where  a  =  27r(5/A  is  a  dimensionless  wavenumber  and  A  is  the  disturbance 
wavelength.  For  reference,  tlu^  Gortler  number  at  the  separation  location  is  also  include<l.  These 
“Gbrtler”  vx)rticos  are  expected  to  be  centered  about  the  outc'i*  region  of  the  velocity  profile,  or 
along 1/2  (Floryan  1986), 

Cro,ssflow  PIV  is  u,sed  to  detect  the  presence  of  streamwise  vortices  by  illuminating  a  cros,s- 
flow  ])lane  normal  to  the  trailing  eclge  surface  13  mm  downstream  from  the  blowing  slot,  or  rOjc 
=  0.025.  Initially,  measurements  are  attempted  at  the  same  test  conditions  used  to  (‘v^luate 
the  Gbrtler  number  in  Figure  40,  Howev^er,  to  boost  the  jet  seed  density  and  image  quality,  the 
chord  Reynolds  number  is  reduced  to  Rcc  =  5.6  x  10'*^  while  maintaining  —  0,014  and  h/c 
=  0.0019.  The  drop  in  Reynolds  number  reduces  the  mass  flow  rate  required  to  maintain  the 
desired  and  hence  increases  the  seed  density,  but  since  the  momentum  coefficient  remains 
constant,  the  same  jet-to- freest  ream  velocity  ratio  is  maintained  =  2).  Counter¬ 

rotating  pairs*  of  vortices,  much  like  those  described  by  Neuendorf  et  al.  (2004),  are  readily* 
visible  in  the  acquired  image  pairs  and  typicalh^  located  between  oiie-half  and  one-and-a-half 
slot  heights  from  the  surface.  An  instantaneous  snapshot  of  the  spanwise  vorticity  computed 
from  an  image  pair  is  presented  in  Figure  41  and  shows  multiple  regions  of  positive  and  negative 


22 


vorticity,  often  located  in  pairs.  A  smaller  region  of  the  flow  is  extracted  in  Figure  42.  wlu‘ie  a 
counter-rotating  pair  of  vortices  is  clearly  present.  The  axes  of  the  vortices  are  located  between 
0.9  <  xjfh  <  1.3,  which,  based  on  the  results  of  the  curved  wall  jet  analysis  for  the  same  and 
=  6.5  X  10’'^,  likely  place  the  vortices  in  the  outer  region  of  the  strearnwise  flow  field. 

Neuendorf  et  aL  (2004)  and  Han  et  aL  (2006)  hypothesize  that  these  vortices  lead  to  separa¬ 
tion  by  “lifting”  the  flow  off  the  surface.  Consider  the  counter-rotating  pair  of  vortices  in  Figure 
42.  Since  these  vortices  are  located  in  the  outer  region  of  U,  they  “pull”  lower- moment  urn  fluid 
from  the  freestream  towards  the  surface  and  “push”  the  higher-moment  inn  fluid  located  near 
yrnax  away  from  the  surface.  In  contrast,  vortex  generators  are  commonly  placed  in  boundary 
layer  flows  to  produce  strearnwise  vorticity  that  “pulls”  the  higher-momentum  freestream  fluid 
towards  the  surface  and  “pushes*’  the  lower- momentum  near-boundary  fluid  away,  delaying  sei>- 
aration.  The  strearnwise  vortices  produced  in  the  outer  region  of  the  curved  wall  jet  do  just  the 
opposite  and  actually  have  implications  to  possible  noise-abating  treatments.  Recall  Slomski 
(2009)  numerically  evaluated  the  effect  of  serrating  the  slot  lip  to  reduce  lip  tones  by  breaking 
up  the  coherent  spanwise  vortices  shed  from  the  lip.  It  is  pos,sible  that  the  saw-tooth  serration 
pattern  could  introduce  these  vortices  upstream  from  where  they  naturally  occur,  causing  the 
flow  to  separate  earlier.  Since  the  computational  domain  evaluated  by  Slomski  did  not  include 
the  separation  region,  the  influence  of  the  serrations  on  flow  separation  could  not  be  deliiiitively 
assess('d. 

4.3  Summary 

Fluid  dynamic  measurements  reveal  a  great  extent  about  both  the  freestream  and  circulation 
control  curved  wall  jet  flows.  Cp  measurements  of  the  circulation  control  airfoil  in  the  open  jet 
test  section  of  the  UFAFF  lack  a  leading  edge  suction  peak  typically  observed  in  closed  test 
section  measurements.  Enclosing  the  test  section  causes  the  leading  edge  peak  to  emerge  but 
does  not  change  the  free.stream  boundary  layer  flow  passing  over  the  slot  lip.  The  influence^ 
of  different  wall  boundary  conditions  on  the  Cp  distribution  is  also  examined  using  data  and 
potential  flow  theory.  The  curved  wall  jet  flow  is  measured  using  PIV  and  found  to  exhibit  sim¬ 
ilarity  in  the  outer  region  of  only.  However,  this  alone  permits  the  development  of  eejuations 
to  predict  the  length  and  velocity  scales  required  for  similarity  as  a  function  of  /?c^.  and 
h/c.  Finally,  flow  separation  is  found  to  be  dependent  on  the  product  CpRcc^  and  strearnwise 
vortices  theorized  to  promote  separation  are  observed  in  crossflow  PIV  meeLsurements. 

The  fluid  dynamic  measurements  also  provide  insight  with  regards  to  acoustics.  Since  the 
turbulent  boundary  layer  passing  over  the  slot  lip  is  independent  of  the  test  section  configuration, 
the  sound  produced  by  the  interaction  of  the  flow  with  the  airfoil  trailing  edge  should  be 
independent  of  the  test  section  configuration  as  well.  The  similarity  of  the  curved  wall  jet  flow 
makes  prediction  of  some  of  the  scales  required  for  How*e’s  model  possible,  particularly  in  the 
assessment  of  curvature  noise  (Howe  2002).  Finally,  serrating  the  slot  lip  for  noisc'  reduction 
may  cause  earlier  flow  separation  by  instigating  the  production  of  strearnwise  vortices  in  the 
outer  region  of  the  flow. 
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5  Acoustics 


The  focus  of  this  section  is  findings  from  acoustics  measurements.  First,  the  presence  and  scaling 
of  tones  is  discussed.  Then,  results  from  phased  array  measurements  arc  presented  in  order  to 
identify  the  primary  noise  sources  under  different  test  conditions.  Different  measurement  and 
processing  techniques  are  used  to  estimate  the  broadband  noise  spectrum.  Finally.  Howe’s 
model  of  circulation  control  acoustics  is  compared  with  a  measurement  (Howe  2002). 

5.1  Tones 

Acoustic  tones  are  undesirable  for  underwater  vehicle  applications  and  a  potential  hurdle  for 
the  application  of  circulation  control  to  underwater  vehicles.  The  presence  of  tones  is  evaluated 
using  microphone  ‘‘Ml”  from  the  test  setup  illustrated  in  Figure  4.  Figure  43  includes  spectra 
for  a  variety  of  momentum  coefficients  at  Rcc  —  6.5  x  lO'^  and  fi/c  =  0.0019.  The  spectra  reveal 
both  low  and  high  frequency  tones  under  certain  conditions.  These  tones  are  e\^luated  further 
in  the  following  sections. 


5.1.1  Low  Frequency  Tones 


Spectra  from  all  test  cases  where  a  low  frequency  tone  is  measured  are  plotted  together  in  F’igurf' 
44.  The  tests  included  cover  variations  in  Rec^  and  h/c.  The  frequency  axis  is  normalized 
using  the  Strouhal  number, 


St 


"^Vte 


(59) 


which  has  a  nominal  value  of  0.21  for  vortex  shedding  from  a  blunt  edge.  The  length  scale  yrt: 
=  4.6  cm,  the  thickness  of  the  trailing  edge  at  the  slot  exit  plane,  is  used  as  the  reference  length 
scale.  The  tones  plotted  in  Figure  44  collapse  at  a  Strouhal  number  just  above  0.21,  indicating 
that  these  low  frequency  tones  are  produced  by  vortex  shedding  from  the  round  airfoil  trailing 
edge. 

With  sufficient  blowing,  these  tones  are  not  measured.  Using  just  a  single  blowing  slot,  the 
tones  are  eliminated  if  >  0.002. 


5.1.2  High  Frequency  Tones 

Previous  research  cited  the  presence  of  high  frequency  tones  thought  to  be  attributed  to  vortex 
shedding  from  the  slot  lip  (Slomski  2009).  Microphones  placed  in  the  UFAFF  measure  high 
frequency  tones  at  a  variety  of  frequencies  and  amplitudes.  A  sampling  of  these  tones  is  provided 
in  Figure  45  for  =  0.014  and  h/c  =  0.0029.  The  frequency  axis  is  again  normalized  by  the 
Strouhal  number,  this  time  defined  by 

(60) 

*  jet 

where  I  is  the  lip  thickness.  The  Strouhal  number  of  the  tones  are  in  close  agreement  near  0.21, 
and  it  is  concluded  that  these  tones  arc  produced  by  vortex  shedding  from  the  slot  lip. 
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Tliese  "lip’'  tones  are  only,  but  non  necessarily,  produced  when  both  freestreain  and  j('t 
Hows  exist.  However,  exact  conditions  corresponding  to  the  generation  of  lip  tone's  could  not 
be  deternimed,  Xo  clear  trends  between  the  appearance  or  niagihtude  of  the  tones  and  How 
parameters  such  as  7?Cc,  h/c,  and  Rcjei  arc  found.  The  directionality  of  the  tones  is  assessed 
to  some  extent  using  the  microphones  placed  on  the  opposite  side  from  the  blowing  slot.  The 
levels  of  the  tones  measured  by  these  microphones  are  at  most  only  slightly  above  broadband 
noise  levels,  as  shown  in  F'igiire  46, 

With  both  vortex  shedding  tones  experimentally  verified,  focus  shifts  towards  broadband 
noise  and  noise  source  identification. 

5.2  Assessment  of  Noise  Sources 

As  discussed  in  Section  3.4,  a  single  microphone  is  unable  to  di.scerii  between  multiple  acoustic 
sources.  Furthermore,  pairs  or  triads  of  microphones  are  unable  to  distinguish  between  multiple 
correlated  noise  sources.  Phased  acoustic  arrays  are,  on  the  other  hand,  distinctly  useful  for 
identifying  noise  sources.  Results  from  phased  array  measurements  of  the  circulation  control 
airfoil  are  presented  in  this  section.  Cases  with  only  freestreain  flow,  only  jet  flow,  and  both 
flows  are  examined. 

5.2.1  In  the  Absence  of  a  Freestream 

The  circulation  control  jot  is  first  considered  in  the  absence  of  a  freestream.  Recall  Figure  !(>, 
which  reveals  that  the  Cp  distributions  for  upper  and  lower  slot  blowing  are  nearly  identical 
in  the  tunnel  open  test  .section  configuration.  It  is  reasonable  to  a.ssiime  that,  since  the  Cp 
distributions  are  similar,  the  trailing  edge  blowing  slot  flows  and  the  sound  they  produce  are 
also  comparable.  Thus,  before  beam  maps  are  discussed,  consider  the  spectra  shown  in  Figure 
47,  presented  as  sound  pressure  level  (SPL)  in  dB,  where 

SPL  =  10 log, 0  •  (fil) 

A/  is  the  bin  width  of  the  processed  spectrum.  These  spectra  are  measured  by  the  center 
microphone  placed  above  the  airfoil  for  cases  where  the  blowing  slot  on  the  s<une  side  as  the 
microphone  is  used,  and  this  is  referred  to  from  this  point  forwards  as  same  side  blowing  (SSB). 
Furthermore,  the  cases  presented  in  Figure  47  are  similar  to  SSB  cases  analyzed  with  the  phased 
array.  Upper  microphone  spectra  are  considered  in  place  of  the  spectra  from  the  array  center 
microphone  since  the  array  microphone  suffers  from  significant  inter-microphone  scattering  at 
higher  frequencies.  The  vertical  lines  in  Figure  47  signify  the  octave  frequencies  where  beam 
maps  are  generated,  notably  1  kHz  to  64  kHz.  Notice  that  when  Rejet  =  1380,  the  sound 
generated  is  hardly  distinguishable  from  the  noise  floor.  Hot-wire  data  used  to  charaderize  slot 
flow  uniformity  indicate  that  the  jet  is  just  becoming  turbulent  at  Rej\>t  =  hUjei/iy  —  2600. 
According  to  Howe  (2002),  slot-jet  interaction  noise  is  eliminated  if  the  jet  is  laminar. 

Beam  maps,  like  the  example  displayed  in  Figure  48,  are  presented  for  a  variety  of  jet 
Reynolds  numbers  in  Figures  49  through  55.  Each  beam  map  includes  overlays  with  the  airfoil. 
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sidewalls,  inlet,  and  diffnser,  as  illustrated  in  Figure  48.  All  nieasurenients  represent  cases  where 
the  blowing  slot  on  the  same  side  of  the  airfoil  as  the  array  is  used.  The  array  point  spread 
function  (PSF)  is  also  shown  in  each  figure  as  reference  for  the  array  resolution.  CSM  diagonal 
removal  is  applied  during  bcamforming  to  eliminate  microphone  self-noise  and  channel  noise 
(Dougherty  2002).  Propagation  times  are  not  corrected  for  shear  layer  refraction  because  of 
the  unknown  effect  of  the  blowing  jet  on  the  tunnel  shear  layer  (Amiet  1978).  Beam  rnap.^  are 
presented  in  absolute  dB  (ref.  20^Pa)  for  comparison  between  different  test  conditions.  The 
scanning  plane  is  located  along  the  chordline,  1.12  in  above  the  array. 

By  and  large,  in  the  absence  of  a  freest  ream,  the  trailing  edge  is  the  dominant  noise  source  at 
all  frequencies,  provided  the  jet  Reynolds  number  is  sufficiently  large.  When  Rej^t  =  1340,  the 
trailing  edge  noise  is  not  uniform  or  no  single  source  is  apparent.  However,  the  levels  are  also 
extremely  low  in  comparison  with  the  higher  cases,  agreeing  with  the  hot-wire  data  that 

indicates  the  jet  is  just  becoming  turbulent  when  Re^et  =  2600.  When  Rej^^i  =  2680,  the  trailing 
edge  is  the  lone  dominant  source,  although  the  sound  power  is  not  distributed  uniformly  along 
the  trailing  edge  at  frequencies  of  4  kHz  and  8  kHz.  In  addition,  the  levels  are  also  extremely 
low  above  16  kHz.  WTien  Re  jet  =  3980  and  higher,  the  trailing  edge  is  the  principal  noise 
source,  the  power  is  generally  distributed  uniformly  along  the  trailing  edge,  and  levels  increase 
with  Re  jet. 

Having  identified  the  trailing  edge  as  the  primary  noi.se  source  in  the  absence  of  a  freestreani, 
the  array  is  used  to  identify  noise  sources  with  the  addition  of  the  tunnel  free.streain. 

5.2.2  With  a  Freestream 

The  addition  of  a  freestream  may  produce  contaminating  noise  sources.  Potentially  harmful 
noise  sources  include  sidewall  scrubbing  noise  and  diffuser  flow'  impingement  noise,  which  was 
severe  enough  in  prior  experiments  in  the  UFAFF,  that  it  motivated  the  remo\^l  of  a  meter-long 
section  of  the  diffuser  to  further  separate  it  from  microphones  and  arrays  (Bahr  2010).  If  these 
noise  sources  are  present,  then  the  array  provides  insight  into  their  strengths  relative  to  the 
sources  of  interest  and  the  frequencies  at  which  they  may  be  problematic. 

Following  the  pattern  of  the  previous  section,  before  beam  maps  are  provided,  consider  the 
upper  microphone  spectra  measured  for  upper  slot  blowing  test  conditions  plotted  in  Figure 
56  for  a  variety  of  inomentuin  coefficients  at  Rcc  —  6.5  x  10*"^.  Like  the  jet-only  spectra,  the 
spectrum  for  Rcjet  =  0  (C^  0)  is  nearly  identical  to  the  Rcjet  =  1330  {Cft  =  0.004)  .spectruni. 

When  Oft  =  0.015  {Rejet  =  2660),  there  is  some  increase  in  SPL  at  frequencies  below  1  kHz 
and  above  10  kHz,  but  in  between  the  difference  is  marginal.  As  increases  further,  the  sound 
levels  increase  over  all  frequencies. 

Beam  maps  are  provided  in  Figures  57  to  63  for  octaves  between  1  kHz  and  64  kHz.  The 
cases  again  represent  SSB.  Diagonal  removal  is  used,  but  a  shear  layer  correction  is  not.  The 
scanning  region  is  located  1.12  in  above  the  array,  along  the  airfoil  chordline.  Array  PSFs  are 
also  shown  for  reference  with  each  figure. 

In  Hgiire  57  (992  Hz),  it  is  immediately  obvious  that  circulation  control  is  not  a  significant 
noise  source  when  the  slot  jet  is  not  turbulent.  Maps  B  through  D  indicate  that  noise  is  coming 
from  the  vicinity  of  the  airfoil,  but  it  is  difficult  to  discern  exactly  where  the  noise  is  originating. 
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There  is  evidence  of  flow  iiiipiiigeiiieiit  noise  on  the  tunnel  diffuser  or  neighboring  ceiling  wedges 
in  map  D,  which  corresponds  to  =  0.017.  At  the  higlic'st  nioineiituin  coefficients  tested,  there 
is  no  evidence  of  flow  iiiipiiigeinent  noise.  In  maps  E  through  G,  the  trailing  edge  region  is  the 
definitive  dominant  noise  source,  although  the  power  appears  skew'ed  to  one  side  of  th(*  trailing 
edge.  The  power  levels  increase  .significantly  with  increasing  beginning  with  —  0.017. 

Maps  produced  at  2  kHz  are  presented  in  Figure  58.  In  maps  B  through  D,  sidewall  scrubbing 
noise,  caused  by  the  turbulent  wall  boundary  layer,  appears  to  be  significant,  and  the  levels 
for  these  maps  are  nearly  identical,  despite  the  fact  that  is  increasing  considerably.  Flow 
impingement  noise  is  also  evident  in  maps  D  and  E  and  is  dominant  in  map  D  where  = 
0.017.  Beginning  with  map  E  {C^i  —  0.037),  the  primary  noise  source  is  located  in  the  vicinity 
of  the  trailing  edge,  but  there  are  large  deviations  in  sound  power  across  the  span  of  the  model 
in  ma{)s  E  through  G. 

The  beam  maps  at  4  kHz  shown  in  Figure  59  have  similar  characteristics  to  the  2  kHz 
maps.  Sidewall  noise  dominates  in  maps  B  through  D  and  is  also  readily  apparent  in  map  E. 
Flow  impingement  noise  is  observed  in  maps  D  and  E  as  well.  In  maps  F  and  G,  the  strongest 
source  appears  along  the  trailing  edge  at  37%  span,  and  its  level  increases  with  C^.  Other,  less 
dominant  sources  are  visible  at  the  junctions  between  the  sidewalls  and  model  trailing  edge. 

In  Figure  60,  which  shows  maps  at  8  kHz,  maps  B  through  D  again  reveal  sidewall  noise,  and 
the  primary  sources  are  located  at  the  intersections  of  the  sidewalls  and  the  airfoil  leading  edge. 
This  sound  could  be  caused  by  horseshoe  vortices  formed  when  the  tunnel  wall  boundary  layer 
passes  around  the  leading  edge.  Compared  with  the  higher  cases,  their  levels  are  relatively 
low.  Ill  inajis  E  through  G,  the  primary  sources  appear  at  the  trailing  edge-sidewall  interfaces 
and  once  again  along  the  trailing  edge  at  37%  sj)an. 

Maps  created  by  the  inner  array  at  16  kHz  are  shown  in  Figure  61.  Sources  in  maps  B 
through  D  include  the  sidewalls,  leading  edge  horseshoe  vortices,  and  diffuser  flow  impingement. 
How'ever,  the  levels  of  these  sources  are  very  low  relative  to  the  sources  at  momentum  coefficients 
of  =  0.037  and  higher.  The  trailing  edge-sidewall  interface  still  dominates  in  map  E.  but 
the  sound  appears  to  distribute  more  uniformly  over  the  trailing  edge  as  increases  in  maps 
F  and  G. 

Beam  maps  for  32  kHz  are  presented  in  Figure  62.  With  no  or  minimal  blowing  (maps 
B  and  C),  the  leading  edge  horseshoe  vortices  and  the  boundary  layer  trip  appear  as  som*ces, 
although  their  levels  are  extremely  low.  The  dominant  sources  are  the  trailing  edge-sidewall 
interfaces  in  map  D  through  G,  and  the  37%  span  location  in  maps  E  through  G. 

Finally,  at  64  kHz  (Figure  63),  maps  A  through  D  are  extremely  noisy,  but  spectra  measured 
by  the  center  microphone  in  Figure  56  indicate  the  levels  at  this  frequency  are  at  or  near  the 
noi.se  floor.  Among  the  sidelohes  apparent  in  maps  E  through  G,  the  dominant  sources  still 
appear  to  be  located  at  the  junctions  of  the  sidewalls  and  the  trailing  edge  and  along  the  trailing 
edge  at  37%  span. 

Figures  57  through  63  reveal  an  assortment  of  noise  sources,  most  of  which  are  uiidesired. 
It  is  clear  that  for  Rec  —  6.5  x  10^,  contaminating  noise  sources  dominate  at  all  frequeiicif»s 
considered  when  <  0.017.  These  sources  include  sidewall  noise,  flow  impingement  noise, 
and  leading  edge  noise,  possibly  from  horseshoe  vortices.  At  =  0.037,  some  of  these  sources 
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are  still  present  with  the  addition  of  noise  at  the  trailing  edge  that  appears  to  originate  at  the 
trailing  edge-sidewall  junctions  and  at  37%  span.  These  same  trailing  edge  sources  appear  to 
dominate  for  the  higher  monientnin  coefficients  considered.  Near  the  sidewall,  a  few  inechaiiisiii.s 
could  be  generating  sound.  First,  vortices  similar  to  wing  tip  vortices  have  been  observed  in 
prior  circulation  control  experiments  with  sidewalls,  and  these  very  vortices  were  the  subject  of  a 
recent  numerical  investigation  (Englar  &  Williams  1972:  Nishino  i:  Shariff  2010).  Second,  there 
is  a  di.scoiitinnity  at  the  edge  of  the  slot  formed  by  the  steel  airfoil  end  plate  just  a  few  millinu'ters 
from  the  foam  wall.  This  edge  could  radiate  sound.  While  these  are  potential  mechanisms  for 
sound  generated  near  the  sidewall,  the  source  appearing  at  37%  span  is  surprising,  especially 
considering  that,  in  the  absence  of  a  freestream,  the  source  is  uniformly  distributed  across  the 
entire  trailing  edge.  Perhaps  the  sound  is  only  generated  with  a  freestream  because  it  is  caused 
by  the  turbulent  boundary  layer  passing  over  the  model  surface  or  lip  edge.  However,  while 
there  is  a  seam  of  two  lip  pieces  near  37%  span,  inspection  reveals  no  distinct  diseontinuity  or 
non-uniformity  eompared  to  the  other  lip  seams,  and  prior  to  the  tunnel  entry,  this  seam  is 
sealed  and  checked  for  leaks. 

To  gather  more  insight  regarding  these  trailing  edge  sources,  beam  maps  from  opposite  side 
blowing  measurements  (OSB)  are  compared  with  SSB  measurements  in  Figures  64  to  70  for 
the  thret^  highest  monientnin  coefficients  considered.  The  left  coliinm  of  beam  maps  (A,  C,  E) 
corresponds  to  SSB,  and  the  right  column  of  beam  maps  (B,  D,  F)  corresponds  to  OSB.  To 
reiterate,  the  array  is  located  on  the  same  side  as  the  blowing  slot  for  the  SSB  cases. 

The  same  and  opposite  slot  blowing  beam  maps  at  992  Hz,  shown  in  Figure  64,  arc  nearly 
undistinguishable,  but  it  is  interesting  that  the  levels  are  slightly  higher  for  OSB.  At  2  kHz 
(Figure  65),  there  are  slight  differences  in  the  sources  identified.  With  OSB,  the  array  does  not 
measure  flow"  impingement  noise  when  =  0.040,  likely  because  the  flow"  impingement  noise 
is  below"  the  array  plane.  At  the  higher  momentum  coefficients,  the  maps  are  fairly  similar,  but 
the  levels  are  as  much  as  5  dB  higher  for  OSB  compared  to  SSB.  At  4  kHz  (Figure  66),  there  are 
substantial  differences  in  the  apparent  noise  sources.  With  OSB,  the  noise  is  much  more  evenly 
distributed  than  SSB,  and  the  levels  are  between  3-4  dB  higher  for  OSB.  At  8  kHz  (Figure  67), 
how"ever,  the  trailing  edge-sidewall  interfaces  begin  to  dominate,  along  with  a  region  near  25% 
span.  Levels  are  just  slightly  higher  for  the  OSB  cases.  It  is  difficult  to  distinguish  the  SSB 
and  OSB  beam  maps  shown  in  F'igure  68  for  16  kHz,  and  for  all  but  the  highest  shown, 
the  maximum  levels  for  SSB  axe  higher  than  OSB.  At  32  kHz  and  64  kHz,  presented  in  Figures 
69  and  70,  respectively,  the  beam  maps  for  OSB  reveal  just  one  source  located  at  one  of  the 
sidewall-trailing  edge  junctions.  The  highest  levels  for  OSB  are  again  3-4  dB  higher  than  SSB 
for  32  kHz,  but  the  maximum  levels  are  higher  for  SSB  at  64  kHz.  It  is  interesting  that  sources 
are  not  observed  along  the  trailing  edge  at  37%  span  and  at  the  opposite  sidew’all  for  OSB. 

The  array  measurements  provide  useful  information  about  the  true  sound  sources  in  this 
experiment.  At  low  momentum  coefficients,  tunnel  sidewall  noise,  flow-  impingement  noise,  lead¬ 
ing  edge  horseshoe  vortices,  and  sidewall-trailing  edge  junction  noise  dominate.  At  the  high^^st 
momentum  coefficients,  sidewall  scrubbing  noise  and  flow  iin{)iiigeuient  noise  are  insignificant, 
but  the  sound  may  largely  be  produced  at  the  sidewall-trailing  edge  interfaces  and  along  the 
trailing  edge  at  37%  span  by  an  unknow’ii  source.  Undoubtedly,  a  single  microphone  is  ctTtainly 
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not  measuring  only  two-dimensional  circulation  control  noise  sources.  Furthermore,  since  the 
trailing  edge  alone  is  never  the  dominant  source,  the  applicability  of  the  coherent  output  power 
and  three- microphone  methods  is  also  questionable. 

5.3  Broadband  Noise 

The  beaiiiforming  results  show  that  multiple  broadband  acoustic  sources  exist  over  all  fre¬ 
quencies  regardless  of  the  iiioineiituni  coefficient  tested.  Since  sound  produccxl  by  the  two- 
dimensional  circulation  control  sources  of  interest  never  appears  to  dominate,  the  COP  and 
three-microphone  methods  iire  likely  not  suitable  for  application  towards  this  investigation. 
However,  they  will  still  be  assessed  and  compared  with  a  single  microphone  autospectnim. 

5.3.1  Free-Standing  Microphones 

The  different  free-standing  microphone  processing  techniques  described  in  Section  3.4.1  are 
compared  for  Rcc  =  6.5  x  10*^,  h/c  =  0.0019,  and  ^  0.10,  which  represents  the  highest 

tested  and,  as  indicated  by  the  beam  maps,  yields  insignificant  sidc'wall  scrubbing  and 
flow  impingement  noise.  Data  are  presented  for  SSB  and  OSB  cases.  The  single  microphone 
spectrum  represents  the  autospectrum  of  the  center  G.R.A.S.  microphone.  The  COP  de- noised 
autospectrum  estimate  is  computed  using  the  center  G.R.A.S.  and  the  downstream  B&K  micro¬ 
phone,  although  the  resultant  spectrum  is  similar  if  the  upstream  microphone  is  used.  Both  the 
COP  and  three- microphone  spectra  represent  estimates  for  the  signal  measured  by  the  center 
G.R.A.S.  microphone. 

Figure  71  compares  the  three  methods  for  SSB.  The  COP  and  thrce-niicrophone  spectra  are 
plotted  as  a  filled  region  that  includes  their  uncertainty  bounds.  The  COP  speetrinn  is  lower 
in  magnitude  than  the  single  and  three- microphone  spectra,  particularly  at  higher  frequencies. 
The  COP  method  is  known  to  under-prediet  levels  if  the  signal-to-noise  ratio  of  the  additional 
microphone  is  not  sufficiently  high  (Bahr  et  al.  2008).  The  thrce-niicrophone  spectrum  is  in 
close  agreement  with  the  single  microphone  spectrum  at  higher  frequencies. 

Similar  results  are  shown  in  Figure  72  for  opposite  side  blowing.  The  COP  spectrum  levels 
are  low  in  comparison  with  the  other  methods,  and  the  three-mierophone  spectrum  agrees  well 
with  the  single  microphone  spectrum,  even  at  higher  frequencies. 

The  general  agreement  of  the  spectra  in  Figures  71  and  72  indicates  that  there  is  little 
uncorrelated  noise  in  the  measurements.  Since  beamforming  results  indicate  the  existence*  of 
multiple  dominant  sources,  these  methods  are  not  appropriate  for  measuring  the  noise  of  interest 
produced  by  the  circulation  control  trailing  edge.  The  three-microphone  and  single-microphone 
spectra  are  in  good  agreement,  suggesting  that  the  thrce-niicrophone  method  provides  a  sat¬ 
isfactory  roprc*sentatioii  of  the  overall  acoustic  field  comprised  of  circulation  control  noise  and 
installation  effects.  Using  the  phased  array,  it  may  be  po.ssible  to  extract  a  spectrum  of  this 
desired  circulation  control  noise  only. 
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5.3.2  Array 


An  alternative  approach  to  the  methods  discussed  in  the  previous  section  is  to  obtain  a  spectrum 
from  the  array  via  integration,  as  described  in  Section  3.4.2.  This  proce.ss  is  certainly  not 
without  its  limitations,  particularly  at  higher  frequencies,  since  the  integrated  levels  are  highly- 
dependent  on  the  accuracy  of  the  array  calibration,  including  the  measured  sensor  locations, 
and  the  region  of  integration. 

An  integrated  spectrum  is  compared  to  the  free-standing  center  G.R.A.S.  microphone  au¬ 
tospectrum  in  Figure  73  for  Rsc  =  6.5  x  10'^,  h/c  =  0.0019,  and  SSB  ^  0.065.  The  outer 
array  is  used  to  compute  the  spectrum  for  frequencies  below  10  kHz,  while  the  inner  array  is 
used  for  frequencies  10  kHz  and  higher.  The  integration  region  is  frequency  dependent  and  aims 
to  exclude  noise  from  the  sidewall-trailing  edge  interfaces.  In  particular,  the  integration  region 
is  defined  by 


BW/2  <  ./■  <  D\V/2 
-L/2  -h  nW/2  <ij<  L/2  -  B\V/2 


where  BW  refers  to  the  computed  array  3  dB  beamwidth,  L  =  1.12  m  is  the  airfoil  si)an,  and 
the  origin  of  the  coordinate  system  is  at  the  midspan  point  where  xfc  =  1.  The  shaded  overlay 
in  Figure  73  represents  the  estimated  uncertainty  in  the  spectral  levels  computed  using  10,000 
iterations  of  the  Monte-Carlo  simulation  described  b}'  Yardibi  et  al.  (2010a).  Both  sj)ectra  in 
Figure  73  are  plotted  using  the  dimensionless  spectrum  form  suggested  by  Howe  (2002), 
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w4iere  U  is  the  freestream  velocity.  4>o  is  the  single-sided  power-spectral  density,  po  is  the 
freest  ream  density,  is  the  source  length,  |x|  is  the  distance  from  the  trailing  edge  at  inidspan 
to  the  observer,  M  is  the  freestream  mach  number,  0  is  the  angle  between  the  chordwise  direction 
and  the  observer,  and  0  is  the  angle  between  the  spaiiwise  direction  and  the  observer.  For  the* 
single  microphone  spectrum,  Lg  =  L  =  l.\2  m,  but  for  the  array  spectrum,  Lg  is  based  on  the 
spanwise  length  of  the  frequency-dependent  integration  region.  As  expected,  the  array  spectrum 
and  single  microphone  spectrum  deviate  significantly,  particularly  at  higher  frequeiich's  where 
the  dominant  sources  appear  at  the  .sidew^all-trailing  edge  interfaces. 

The  integrated  array  spectrum  is  used  in  the  following  section  for  comparison  with  the 
acoustic  model  derived  by  Howe  (2002). 


5.4  Howe’s  Model  of  Circulation  Control  Acoustics 

Howe  (2002)  theorized  three  significant  broadband  noi.se  types:  curvature  noise,  pa.ssiv('  slot 
noise,  and  slot-jet  interaction  noise.  Curvature  noise  is  produced  by  boundary  layer  turbulence 
scattering  off  the  rounded  trailing  edge.  Passive  slot  noise  is  generated  by  freestream  boundary 
layer  turbulence  scattering  off  the  slot  lip.  Finally,  slot-jet  interaction  noise  is  caused  by  the 
interaction  of  turbulence  in  the  trailing  jet  with  the  slot  lip  and  trailing  edge  surface  near  the 
slot.  Howe  derived  mathematical  models  for  these  three  noise  sources  that  are  functions  of 
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local  How  scales,  including  the  displacement  thickness,  iiieaii  velocity,  and  friction  velocity.  In 
order  to  accurately  compare  Howe’s  model  with  the  experimental  results,  these  scales  m'ed  to 
be  determined.  First,  the  primary  assumptions  behind  Howe’s  model  are  introducixi. 

5.4.1  Assumptions 

Howe’s  model  follows  three  primary  assumptions.  First,  the  freestream  Mach  number.  M  1. 
For  the  cases  considered,  M  =  0.06,  so  this  assumption  is  satisfied.  Second,  the  frt'ostream 
Reynolds  number  must  be  large.  Since  Rcc  —  6.5  x  10*^,  this  assumption  is  also  satisfied. 
Finally,  Howe  assumed  compactness  of  the  scattering  edge  thickness.  For  curvature  noi.se,  Howe 
considered  the  trailing  edge  radius  as  the  length  scale  of  interest,  and  hence,  R  X.  For  the 
circulation  control  airfoil  of  the  present  study,  R  =  X  when  /=  15  kHz,  so  provided  /  <  1  kHz, 
the  compactness  assumption  is  met.  Howe  considered  the  lip  thickness  as  the  length  scale  of 
interest  when  determining  compactness  for  passive  slot  noise  and  slot-jet  interaction  noise.  The 
airfoil  under  investigation  has  a  lip  thickness  of  0.28  mm.  An  acoustic  wavelength  of  0.28  inin 
corresponds  to  a  frequency  of  over  1200  kHz,  so  this  compactness  assumption  is  also  valid  for 
the  highest  frequency  of  interest,  80  kHz. 

5.4.2  Estimates  of  Flow  Scales 

Howie’s  model  is  assessed  using  flow  meeusurements  for  Rec  =  6.5  x  10^,  C^i  —  0.057  (Rejt^t  = 
5000),  and  h/c  =  0.0019.  This  case  represents  the  highest  momentum  coefficient  studied  with 
PIV.  The  flow'  scales  (mean  velocity,  displacement  thickness,  and  friction  velocity)  required  to 
evaluate  Howe’s  model  are  estimated  from  this  data. 

Curvature  noise  Howe  (2002)  defines  curvature  noise  as  the  sound  produced  by  turbulence  in 
the  boundary  layer  piissing  over  the  rounded  trailing  edge.  This  is  perhaps  an  unclear  definition, 
since,  as  sliowii  in  Section  4,  the  flow  over  the  rounded  trailing  edge  is  that  of  a  w'all  jet  and 
not  a  conventional  boundary  layer.  In  addition,  exactly  where  to  estimate  the  scales  along  the 
trailing  edge  is  vague.  In  his  analysis,  Howe  used  the  point  of  maximum  eurwature  on  an  ellipse. 
Since  the  airfoil  trailing  edge  under  investigation  has  a  constant  radius,  curvature  is  constant. 
Instead,  since  it  is  well-known  that  turbulent  sound  production  increases  with  velocity,  scales 
for  eurvatiire  noise  are  evaluated  at  the  point  of  maximum  velocity  along  the  trailing  edge. 

Th('  mean  tangential  velocity  |)rofiie  at  the  position  of  maximum  velocity,  found  to  be  r6/c  = 
0.0129.  is  plotted  in  Figure  74.  The  three  scales  of  interest  include  the  mean  velocity  outside  the 
boundary  layer,  the  displacement  thickness,  and  the  friction  velocity.  Since  the  flow  represents  a 
wall  jet,  the  inner  region  provides  the  closest  resemblance  to  a  conventional  boundary  layer  and 
will  henceforth  be  used  to  determine  these  flow  scales.  The  mean  velocity  outside  the  boundary 
layer  is  chosen  to  l)e  the  local  maximum  velocity,  which  is  measured  to  be  Uq  =  U„iaj  =91.3 
m/s.  The  displacement  thickness  and  friction  velocity  are  much  more  difficult  to  calculate 
because  of  insuffieient  data  near  the  surface.  Furthermore,  the  maxiininn  velocity  is  found  to 
be  the  closest  datum  point  to  the  surface.  As  an  approximation,  a  one-seventh  powvr  law  is 
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used  to  estimate  the  displacement  thickness  and  momentum  thickness  (Prandtl  19G1). 
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are  determined  by  evaluating  the  power  law.  The  skin  friction  coefficient  is  then  estimated  using 
the  Karman  integral  relation  (White  2006), 
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where  H  =  5^/0  \s  the  shape  factor,  and  Re^  =  Umax^/^  the  Reynolds  number  based  on  the 
iiiomeiiturn  thickness  and  '^freestreain'  velocity.  After  the  skin  friction  is  obtained,  the  wall 
shear  stress  and  friction  velocity  are  estimated. 

Tw  =  (57) 

(68) 

The  displacement  thickness  and  friction  velocity  are  approximately  0.010  nirii  and  9.57  in/s, 
respectively. 


Passive  slot  noise  Howe  (2002)  considers  passive  slot  noise  to  be  the  sound  produced  by 
turbulence  in  the  exterior  freestream  boundary  layer  scattering  off  the  slot  lip.  Like  curvature 
noise,  the  as.sumption  of  a  conventional  boundary  layer  flow  pas.sing  over  the  lip  is  not  entirely 
true.  Figure  22  shows  that  the  flow  passing  over  the  lip  edge  takes  the  form  of  a  wall  jet 
more  so  than  a  boundary  layer.  Like  curvature  noise,  the  mean  velocity  outside  the  bound<u:y 
layer  is  taken  as  the  local  maximum  velocity,  and  the  displacement  thickness  is  computed  using 
Equation  64. 

The  mean  tangential  velocity  profile  measured  using  PIV  is  shown  in  Figure  75.  The  max- 
iiniim  velocity  is  found  to  be  Us  =  Umax  —  41.4  m/s.  Becau.se  there  is  sufficient  data  in  the 
profile,  the  displacement  and  momentum  thicknesses  are  computed  using  the  data  alone.  The 
friction  velocity  is  estimated  using  the  Karman  integral  relation.  The  displacement  thickness  is 
0.140  mm,  and  the  friction  velocity  is  2.78  m/s. 


Slot-jet  interaction  noise  Slot-jet  interaction  noise  is  theorized  by  Howe  (2002)  to  he  gen¬ 
erated  by  the  interaction  of  turbulence  in  the  jet  with  the  slot  lip  and  lower  slot  surface,  i.e. 
the  trailing  edge  in  the  vicinity  of  the  slot.  To  estimate  the  flow  scales  for  this  noise  source, 
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the  mean  tangential  velocity  profile  just  downstream  of  the  slot  exit  plane  is  evaluated.  The 
mean  velocity  is  once  again  taken  to  be  the  local  maximum  wall  jet  velocity.  Because  the  uj)per 
portion  of  the  profile  is  influenced  by  tlu*  lip  wake,  the  displacreiiieiit  and  moment um  thicknesses 
are  analyzed  at  the  lower  Coanda  surface  using  the  one-seventh  power  law  of  Equation  63.  The 
Karman  integral  relation  is  again  used  to  estimate  the  friction  velocity. 

The  mean  tangential  velocity  profile  at  the  slot  exit  is  shown  in  Figure  76.  The  maximum 
velocity  is  found  to  be  !7j  =  f/mox  =  88.8  m/s,  the  displacement  thickness  is  estimated  as  0.025 
mm,  and  the  friction  velocity  is  approximated  to  be  6.84  m/s. 

5.4.3  Comparison  with  Measurement 

Now  that  scales  have  been  estimated  from  the  PIV  data,  they  will  used  to  be  compare  Howe's 
model  to  a  measured  spectrum.  For  conii)arison  puri)oses,  first,  Howe’s  model  is  assessed  using 
the  scales  it  estimates  given  just  the  model  geometry,  freest  ream  velocity,  and  jet  velocity. 
SpcH'ifically,  the  scales  are  the  freestream  velocity,  U,  the  jet  velocity  Ujet.  the  semi-major  axis 
of  the  airfoil  elliptic  profile,  «,  the  trailing  edge  radius,  r,  the  slot  height,  h,  and  the  lip  thickness, 
1. 


U  =  18.9  m/s,  =  80.6  m/s,  a  =  0.2606  in 

r  =  0.0222  m.  /i  =  1  x  10“^  m,  I  ^  2.1  x  lO""*  m 

The  spectrum  produced  by  Howe’s  model  for  these  inputs  is  shown  in  Figure  77  along  with  the 
spectrum  obtained  from  the  array  measurement.  The  spectra  are  presented  in  the  form  given 
by  Equation  62.  For  the  present  case,  values  for  these  constants  are  given  below. 

A)  =  1.21  kg/m^,  |x|  =  1.12  m,  M  =  U/cq  =  0.055 
0  =  7r/2  rad,  V-’  =  7’’/2  rcicl 

Recall,  when  applied  to  the  array  spectrum,  L*.,  the  source  length,  is  frequency-dependent  and 
based  on  the  array’s  3  dB  beamwidth.  For  the  predicted  spectrum  in  Figure  77,  —  1.12  ni. 

As  evident  in  Figure  77,  there  are  significant  differences  betweem  the  measured  spectra 
and  Howe’s  model.  However,  recall  Howe’s  spectrum  is  computed  using  only  estimates  of 
the  length  and  velocity  scales.  These  scales  are  listwl  in  Table  4  and  compared  with  those 
determined  from  the  PIV  measurements.  Not  surprisingly,  there  are  large  differences  bctwcHUi 
the  scales,  especially  those  for  curvature  noise.  Recall,  Howe  (2002)  considers  curvature  noise  to 
be  produced  by  boundary  layer  turbulence,  whereas  in  the  actual  flow,  there  is  a  wall  jet  instead 
of  a  conventional  boundary  layer  passing  over  the  trailing  edge.  If  the  inner  region  of  the  wall 
jet  is  taken  to  be  the  boundary  layer,  then  the  mean  and  friction  velocities  will  certainly  be 
higher,  and  the  displacement  thickness  will  be  significantly  smaller. 

When  the  scales  estimated  from  the  PIV  measurements  are  substituted  into  Howe's  model, 
the  predicted  spectrum  is  quite  different,  as  shown  in  Figure  78.  The  differences  betwt^'ii  the 
measured  and  predicted  spectra,  though,  are  still  significant.  The  shaded  regions  presented  with 
each  predicted  spectrum  represent  uncertainty  bounds,  estimated  using  10.000  iterations  of  a 
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Monte-Carlo  simulation.  The  friction  velocities,  displacement  thicknesses,  moan  velocities,  slot 
height,  lip  thickness,  and  freestreani  velocity  are  all  perturbed  with  each  iteration.  Tin*  pertur¬ 
bations  to  the  friction  velocities,  displacement  thicknesses,  and  mean  velocities  are  themsehTs 
determined  from  10,000  iterations  of  Monte-Carlo  simulations  based  on  the  PIV  measurement 
uncertainty.  Even  considering  all  uncertainties,  the  differences  between  the  predicted  and  nieci- 
sured  spectra  are  large.  However,  some  useful  information  may  still  be  extracted  from  Howe's 
model.  According  to  the  model,  passive-slot  noise  is  the  dominant  noise  source  for  frequencies 
below  20  kHz,  and  slot-jet  interaction  noise  is  the  primary  source  at  higher  frequencies. 

Finally,  an  attempt  is  made  to  determine  values  for  the  flow  scales  in  order  for  Howe's 
model  spectrum  to  approach  the  measured  array  spectrum.  Howe’s  model  is  fit  to  the  measured 
spectrum  using  a  least  squares  approach,  and  the  flow  scales  that  provide  the  best  fit  are 
determined.  The  results  from  the  least  squares  fit  are  dependent  on  the  initial  guess  and  the 
bounds  set  for  each  parameter.  The  initial  guesses  for  the  parameters  are  the  measured  values  of 
the  scales  listed  in  Table  4.  The  lower  bounds  for  all  parameters  are  set  to  0.  The  upper  bounds 
for  the  mean  velocities  are  set  to  100  m/s,  the  upper  bounds  for  the  friction  velocities  are  set  to 
20  m/s,  and  the  displacement  thickness  bounds  for  the  wall  jet  and  lip  boundary  layer  profiles 
are  set  to  0.5  mm  and  1  mm.  respectively.  The  best  fit  spectrum  is  (oinpared  to  the  measured 
spectrum  in  Figure  79,  where  the  two  are  shown  to  be  in  close  agreement.  Notice  there  is  no 
contribution  from  slot-jet  intcreu^tion  noise.  To  reiterate,  the  resultant  Ih'sI  fit  paranu’ters  are 
dependent  on  the  initial  guesses  and  bounds,  so  this  solution  is  provided  simply  for  illustrative 
purposes.  The  scales  determined  from  the  best  fit  are  listed  in  Table  5.  There  are  generally 
large  variations  between  the  scales  obtained  from  the  PIV  measurements  and  the  scales  listtni 
in  Table  5.  In  addition,  some  of  the  scales  are  unreasonable,  including  u*  =  19  m/s  for  the 
flow  passing  over  the  lip,  =  0.5  mm  with  a  corresponding  i;*  =  20  m/s  for  the  wall  jet  flow 
responsible  for  curvature  noise,  and  =  1.0x10“^  mm  for  the  displacement  thickness  at  the 
slot  exit. 

5.5  Summary 

The  sound  produced  by  a  circulation  control  airfoil  in  an  open  jet  anechoic  wind  tunnel  is 
measured  and  characterized  using  free-standing  microphones  and  a  nest(‘d  phased  acoustic  ar¬ 
ray.  Tones  produced  by  vortex  shedding  from  the  round  trailing  edge  and  blunt  slot  lip  are 
detected.  Beainfonning  maps  indicate  that  the  trailing  edge  is  the  dominant  source  in  the 
absence  of  a  freestream,  but  with  a  freestreani,  there  are  a  multitude  of  sources.  At  low  nio- 
mentuin  coefficients  corresponding  to  laminar  slot  jets,  contaminating  noise  source's,  including 
flow  impingenient  noise  and  sidewall  scrubbing  noise,  are  observed  to  dominate.  At  higher  mo 
inentum  coefficients,  noise  from  the  trailing  edge  is  the  primary  source,  but  the  sound  power 
is  unevenly  distributed  along  the  trailing  edge.  The  beam  maps  indicate  that  the  interfaces 
of  the  sidewalls  and  trailing  edge  are  the  primary  sources,  particularly  at  higher  frequencies. 
Spectra  computed  using  the  COP  and  three-microphone  methods  are  compared  with  single 
microphone  autospectrum  and  found  to  be  in  close  agreement,  indicating  little  imcorrelattHl 
noise  in  the  measurements.  All  three  methods  are,  by  their  nature,  unable  to  distinguish  the 
twodimensional  circulation  control  noise  generated  at  the  trailing  edge  with  noise  produced  by 
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the  aforementioned  contaminating  noise  sources.  A  spectrum  is  also  computed  using  the  phased 
array  by  integrating  over  a  small,  frequency-dependent  region  of  the  trailing  edge  to  ininiinize 
the  sidewall  end  effects.  This  integrated  spectrum  is  compared  with  a  model  of  ciiculation 
control  acoustics  presented  by  Howe  (2002)  and  fotmd  to  differ  on  average  by  30  dB  How('*s 
model  suggests  that,  for  the  test  conditions  considered,  passive  slot  noise  is  the  dominant  noise 
source  for  frequencies  below  20  kHz,  and  slot-jet  interaction  noise  is  the  principal  source  at 
higher  frequencies. 

6  Conclusions  and  Future  Work 

6.1  Key  Findings 

6.1.1  Fluid  Dynamics 

Midspan  surface  pressure  measurements  of  an  elliptic  circulation  control  airfoil  in  the  UFAFF 
reveal  the  absence  of  a  leading  edge  suction  peak  regularly  observed  in  circulation  control 
experiments  and  numerical  studies.  The  lack  of  a  leading  edge  suction  peak  contributes  to  a 
deficit  in  the  lift  produced  compared  with  prior  experiments  of  the  same  ge^Dmotry  in  closed 
test  section  wind  tunnels.  The  leading  edge  suction  peak  is  recovered  by  enclosing  the  wind 
tunnel  test  section,  and  the  lift  deficit  is  eliminated.  However,  the  surface  pressure  distribution 
is  highly-dependent  on  the  boundary  conditions  enforced  by  the  suction-side  tunnel  boundary. 
Measured  surface  pressure  distributions  are  compared  with  potential  flow  theory  for  flow  around 
an  ellipse  using  the  conformal  mapping  technique.  The  theory  provides  insight  into  the  behavior 
of  the  flow  and  explains  the  elevated  snetion-side  and  reduced  pressure-side  surface  pressure 
magnitudes  measured  with  a  rigid  versus  porous  suction-side  tunnel  boundary.  In  addition,  the 
theory  also  supports  the  observations  that  the  pressure-side  tunnel  boundary  is  nearly  irrelevant 
in  determining  the  surface  pressure  distribution. 

PIV  measureiiients  reveal  the  extent  to  which  eiielo.sing  the  test  section  significantly  modific'^ 
the  loading  edge  flow  field.  In  an  open  jet  test  .sec  t ion,  leading  edge  stagnation  point  movememt 
is  minimal  with  increasing  momentum  coefficient.  In  a  closed  t(\st  .sc'ction,  leading  edge  stag¬ 
nation  point  movement  is  considerably  more  significant  as  the  momentum  eoefficicmt  inerea.s(\s. 
Contrary  to  the  leading  edge  observ'ations,  PIV  measurements  reveal  negligible  difference's  in 
the  boundary  layer  flow  passing  over  the  slot  lip  when  the  test  section  is  opened  or  closed.  The 
scattering  of  turbulent  pressure  fluctuations  in  this  boundary  layer  off  the  slot  lip  is  one  of  the 
noise  sources  theorized  by  Howe  (2002).  Hence,  this  noise  source  mechanism  should  remain 
largely  unchanged  regardless  of  whether  the  test  section  is  enclosed  or  not.  This  finding  can 
provide  direction  to  future  researchers  considering  circulation  control  experiments  in  open  or 
closed  wind  tunnels  when  acoustic  measurements  may  be  of  interest. 

Additional  PIV  measurements  focus  on  the  curved  wall  jet  flow  and  its  similarity.  Although 
a  similarity  solution  indicates  that  full  flow  similarity  is  only  achievable  if  the  curved  surface 
takes  the  shape  of  a  logarithmic  spiral,  the  outer  regions  of  the  mean  tangential  velocity  profile.s 
do  indeed  exhibit  similarity  u.sing  scales  based  on  the  maxinium  velocity  and  the  Reynolds  stress 
or  mean  shear  (Launder  &  Rodi  1983).  The  length  and  velocity  scales  required  for  similarity  are 
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ineasiired  for  a  collection  of  test  caries  where  the  chord  Reynolds  nuinber,  nioiiientinn  coeflicieiit, 
and  slot  height  are  varied.  The  data  for  these  scales  is  assembled  and  found  to  collapse,  and 
the  resultant  best-fit  equations  of  the  scales  are  a  function  of  the  product  of  the  chord  Reynolds 
number  and  momentum  coefficient.  This  so-called  Reynolds  corrected  momentum  coefficient 
can  also  be  recast  as  a  product  of  the  jet  Reynolds  number  and  Jet-to-freestrcam  velocity  ratio 
(see  Equation  54)  indicating  that,  unlike  the  case  of  the  planar  wall  jet  in  a  freestreani,  the 
length  and  velocity  scales  describing  the  curved  wall  jet  flow  in  a  freestreain  are  dependent  on 
the  jet  Reynolds  number  (Zhou  &  Wygnanski  1993).  The  equations  for  predicting  these  length 
and  velocity  scales,  given  by  Equations  50  through  53,  compjire  favorably  with  measurements. 

Flow  separation  is  also  asse.ssed  using  PIV,  and,  like  the  length  and  velocity  scales  of  the 
flow,  the  separation  location  is  a  function  of  the  Reynolds  corrected  momentum  coefficient. 
An  equation  for  predicting  the  separation  location  is  provided  and  found  to  match  the  present 
measurements  with  excellent  accuracy.  However,  the  equation  may  only  be  valid  for  the  airfoil 
geometry  investigated,  as  using  it  to  predict  the  separation  location  for  a  different  geometry 
studied  by  Novak  &  Cornelius  (1986)  and  Novak  et  al.  (1987)  leads  to  large  errors.  Regarding 
the  separation  mechanism,  crossflow  PIV  measurements  uncovered  streainvvise  pairs  of  counter¬ 
rotating  vortices  only  previously  measured  in  a  curved  wall  jet  in  the  absence  of  a  freestreain 
(Likhachev  et  al.  2001;  Neuendorf  et  al  2004;  Han  et  al  2006).  These  vortices,  which  “pulF' 
high-monieiit  11111  fluid  away  from  the  surface,  are  thought  to  he  ultimately  responsihle  for  flow 
separation  (Neuendorf  et  al.  2004;  Han  et  al.  2006).  If  that  is  the  case,  then  these  vortices 
may  preclude  the  use  of  serrated  lip  treatments  to  reduce  lip  vortex  shedding  tones,  since  the 
serrations  may  stimulate  and  advance  the  earlier  development  of  these  vortices  in  the  flow. 

6.1.2  Acoustics 

Acoustic  measurements  reveal  the  presence  of  low  and  high  frequency  tones  associated  with 
vortex  shedding  from  the  round  trailing  edge  and  blunt  slot  lip,  respectively.  The  low  frequency 
tone  can  be  eliminated  with  single  slot  blowing  provided  the  momentum  coefficient  is  at  least 
0.002.  The  emergence  and  behavior  of  the  high  frequency  tones  is  somewhat  sporadic,  as  no 
clear  trends  biised  on  flow  parameters  like  the  chord  Reynolds  number,  nionuaiturn  coefficient, 
jet  Reynolds  number,  and  slot  height-to-chord  ratio  are  found.  However,  the  lip  tones  are  only, 
but  not  nece.ssarily,  produced  when  both  the  chord  and  jet  Reynolds  numbers  are  nonzero.  Also, 
the  tones  are  hardly  discernible  from  the  broadband  noise  when  measured  on  the  opposite  side 
from  the  blowing  slot,  indicating  that  these  high  frequency  tones  are  highly  directive. 

Broadband  noise  sources  are  identified  using  a  nested  phased  microphone  array.  Without  a 
freestreain,  the  trailing  edge  is  the  dominant  sound  source  provided  the  jet  Reynolds  number 
is  sufficiently  high  to  create  a  turbulent  jet.  With  a  freest  ream,  however,  many  sources  are 
measured.  At  the  chord  Reynolds  number  tested,  Rec  =  6.5  x  lO'*^,  containinatiiig  noise  from 
sidewall  scrubbing  and  flow  impingement  are  found  to  dominate  at  low  momentum  coefficients  or 
jet  Reynolds  numbers.  In  addition,  leading  edge  horseshoe  vortices  formed  at  the  sidewall- model 
junctions  dominate  at  frequencies  below  8  kHz  for  <  0.017,  although  their  levels  are  low 
relative  to  the  other  sources  identified.  At  higher  momentum  coefficients,  the  trailing  edge  is  the 
dominant  source.  However,  the  sound  is  not  evenly-distributed  along  the  trailing  edge.  Instead. 
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sources  are  measured  at  the  trailing  e^lge-sidewall  junctions,  likely  caused  by  the  interaction 
of  induced  streamwise  vortices  with  the  sidewalls,  and  near  37%  span,  the  origin  of  which  is 
unknown  but  may  be  due  to  the  seam  between  two  airfoil  components.  Array  measurenierits 
of  tests  utilizing  the  opposit(?  blowing  slot  from  the  array  indicate  that  at  frequencies  of  8  kHz 
and  higher,  the  trailing  edge-sidewall  junctions  are  the  primary  sources. 

The  presence  of  multiple  sources  makes  the  estimation  of  broadband  nois(3  produced  by  the 
circulation  control  airfoil  extremely  difficult.  Spectra  computed  using  different  microphone  pro¬ 
cessing  techniques,  including  single  microphone  autospectrum,  coherent  output  power,  tand  the 
three- microphone  method,  are  found  to  be  in  close  agreement.  However,  since  the  array  (iete'<'ts 
multiple  sources  over  all  frequencies  of  interest,  tlie  results  from  the  microphone  processing 
techniques  are  all  invalid  since  the  single-source  assumption  is  violated.  Thus,  these  methods 
are  not  suitable  for  measurements  of  “two-dimensionaf'  circulation  control  noise,  including  the 
sources  theorized  by  Howe  (2002).  Array  integration  techniques  provide  the  only  means  to 
estimate  the  sound  not  attributed  to  the  undesired  noise  sources,  but  more  work  is  neiKled  to 
perfect  calibration  techniques,  and  eissess  and  reduce  microphone  scattering  effects. 

Finally,  Howe's  model  of  circulation  control  acoustics  is  evaluated  at  Rcc  =  6.5  x  I0‘\ 

=  0.057,  and  h/c  =  0.0019  where  circulation  control  related  noise  is  dominant  (Howe  2002). 
The  mean  velocity,  friction  velocity,  and  displacement  thickness  scales  required  by  Howe’s  model 
are  estimated  from  the  PIV  ineasureiiients  of  the  trailing  edge  How  field.  Using  these  scales, 
large  differences  (~  30  dB)  are  observed  between  Howe’s  model  and  a  broadband  spectrum  ol> 
tained  from  array  measurements.  Howe's  model  indicates  that,  at  the  test  conditions  evaluated, 
which  correspond  to  a  low  chord  Reynolds  number  and  a  moderate  momentum  coefficient,  the 
dominant  noise  mechanism  for  frequencies  below  20  kHz  is  the  turbulent  freestreain  boundary 
layer  flow  passing  over  the  slot  lip.  The  interaction  of  turbulence  in  the  slot  jet  with  the  slot 
lip  is  the  dominant  source  at  higher  frequencies.  Noise  produced  by  the  turbulent  flow  passing 
over  the  round  trailing  edge  is  negligible. 

6.2  Research  Impact 

The  initial  goals  of  this  research  investigation  placed  an  emphasis  on  identifying  and  charac¬ 
terizing  the  two-dimensional  noise  source  mechanisms  of  a  circulation  control  airfoil.  However, 
the  rather  unexpected  results  of  this  study  indicate  that  three-dimensional  noise  sources  may 
be  of  more  interest  and,  from  an  acoustics  perspectiv^e,  a  primary  deterrent  to  the  application 
of  circulation  control  to  underwater  v'ehicles.  The  noise  .sources  identified  at  the  junctions  b<v 
tween  the  trailing  edge  and  sidewalls  are  significant  at  frequencies  of  at  least  4  kHz  and  higher. 
These  noise  sources,  not  the  original  two-dimensional  noise  sources  considered,  warrant  further 
investigation. 

The  findings  of  this  iiwestigation  also  highlight  a  recurring  theme  in  experimental  aeroa- 
constics  -  quantifying  the  source  of  interest  is  often  extremely  difficult.  The  noi.se  produced 
by  the  baseline  or  low  slot  blowing  flows  is  heavily  contaminated  by  sidewall  scrubbing  and 
flow  impingeinent  noise.  Those  imdesirc'd  sources  must  In'  eliminated  if  the  noise  produced  by 
circulation  control  is  to  be  as.sessed  at  such  test  conditions.  Of  course,  that  is  not  a  straight¬ 
forward  task,  and  options  are  limited  based  on  the  flow  facility  and  the  model  being  tested. 
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For  example,  if  sidewalls  are  eliminated,  other  uiidesired  sources  or  three-dimeiisioiiHl  effects 
would  be  introduced.  Perhaps  a  sidewall  suction  system,  if  inipleinented  in  a  quiet  fashion, 
could  eliminate  the  three-dimensional  effects  without  introducing  additional  noise.  As  for  nu‘a- 
surement  techniques,  the  multiple  sources  identified  by  the  phased  acoustic  array  prohibit  the 
use  of  multiple  microphone-based  methods  if  the  two-dimensional  sound  sources  are  of  primary 
interest.  However,  these  methods  might  be  suitable  in  the  study  of  the  three-dimensional  noise 
sources.  Otherwise,  array-based  methods  provide  the  only  possible  manner  by  which  to  study 
the  two-dimensional  sources  of  interest. 

6.3  Recommendations  for  Future  Work 

Since  this  investigation  has  identified  substantial  three-dimensional  noise  source  mechanisms, 
future  wind  tunnel  measurements  should  focus  on  three-dimensional  models  similar  to  an  un¬ 
derwater  vehicle  control  surface.  A  scaled  semi-span  circulation  control  wing  should  be  mounted 
to  a  section  of  a  model  veliicle  hull  and  tested  in  an  anechoic  wind  tunnel.  Furthermore,  by 
removing  the  porous  tunnel  walls,  sidewall  scrubbing  noise  will  be  eliminated.  Instead,  the  rela¬ 
tive  strength  of  the  sound  produced  by  the  hull-trailing  edge  noise  source  can  be  evjiluated,  aiul 
other  three-dimensional  effects  representative  of  an  actual  underwater  vehicle  can  he  gauged. 
Ill  addition,  future  models  should  be  constructed  in  as  few  components  as  j)ossible,  eliminating 
the  likelihood  that  model  seams  or  joints  can  influence  the  measurements. 

Additional  work  is  needed  to  fully  characterize  the  emergence  and  behavior  of  the  high- 
frequency  lip  tones.  These  tones  could  be  a  significant  impediment  to  the  application  of  cir¬ 
culation  control  to  underwater  vehicles.  Since  the  tones  have  been  experimentally  verified,  lip 
modifications  aimed  at  reducing  their  levels,  like  those  suggested  by  SloiiLski  (2009),  should  be 
experimentally  evaluated.  Furthermore,  the  effect  of  such  modifications  on  flow  separation  must 
be  adequately  assessed. 

The  present  investigation  utilizes  PIV  to  measure  the  length  and  velocity  scales  re'quired 
for  Howe’s  model,  but  PIV  is  unable  to  measure  extremely  close  to  surfaces,  which  is  t^ssential 
for  highly-accurate  displacement  thickness  and  friction  velocity  estimates  of  the  curved  wall 
jet  (Howe  2002).  Additional  measurements  are  recommended  on  large r-scale  models  in  larger 
test  facilities,  so  that  spatial  resolution  is  subsequently  improved.  Hot-wire  anemometry  should 
also  be  considered  in  place  of  PIV  for  measuring  near- wall  velocity  profiles,  provided  the  flow* 
is  nominally  unidirectional.  Microelectromechanical  systems  (MEMS)  shear  stress  sensors,  like 
the  one  developed  by  Chandrasekhar  an  (2009),  should  be  considered  for  installation  along  the 
Coanda  surface.  Although  packaging  such  devices  may  be  difficult  on  a  highly-contoured  .sur¬ 
face,  they  could  provide  a  precise  measurement  of  the  friction  velocity.  Alternatively,  global 
interferometric  skin-friction  measurements  of  the  Coanda  surface  and  lip  could  also  provide 
accurate  friction  velocity  information  (Xaughton  k  Sheplak  2002).  Highly-accurate  flow^  mea¬ 
surements  at  a  variety  of  test  conditions  are  essential  to  confirm  under  w'hat  conditions  passive 
slot  noise  is  indeed  the  primary'  two-dimensional  noise  source. 

Finally,  further  work  is  needed  to  support  the  use  of  array-based  methods  in  aeroacoustic 
measurements.  These  methods  are  particularly  challenging  when,  like  in  the  present  investi¬ 
gation,  high  frequency  sources  must  be  considered.  Higher  frequencies  motivate  the  need  for 
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smaller  aperture  arrays,  but  smaller  aperture  arrays  introduce  significant  iniciophoiie  scattering 
effects.  Either  the  scattering  effects  nc'ed  to  he  characterized,  or  a  suitable  treatnient  needs  to 
be  found  to  mitigate  scattering.  The  latter  may  not  be  too  effective  in  some  circumstances, 
like  the  present  investigation,  where  the  distance  separating  microphones  is  less  than  a  half¬ 
wavelength  at  the  highest  frequency  of  interest.  Instead,  an  array  of  MEMS  microphones  should 
be  considered  for  future  experiments.  Such  an  array  would  permit  closer  sensor  spacing  and  pro¬ 
vide  reduced  scattering  effects,  both  from  neighboring  sensors  and  the  microphones  themselves 
(Underbrink  2002;  Arnold  et  ai  2003). 
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Figure  1:  Circulation  control  airfoil  geometry,  a)  Overall  dimensions,  b)  Trailing  edge  gtx)iiietry 
(/  is  the  lip  thickness,  and  t  is  the  ellipse  thickness). 
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Figure  2:  Circulation  control  airfoil  installed  in  the  UFAFF.  a)  2009.  V))  2010. 
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FiRuro  3:  PIV  trfiiling  edge  moa,suromont  regions:  (a)  trailing  edge  wake,  (b)  flow  over  lip,  (f 
curved  wall  jet,  and  (d)  curved  wall  jet  separation.  Airfoil  geometry  sketched  for  reference. 
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Figure  4:  Microphone  experimental  setup. 
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Figure  5:  Phased  array  experimental  setup. 


Figure  6;  Scanning  region  of  I  =  1 . L  points. 
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(a)  (c) 


Figure  7:  Pictures  of  the  phased  array,  a)  Installed  in  the  UFAFF,  looking  towards  the  inlet, 
b)  Magnified  view  of  the  array,  c)  Acoustic  treatment  removed  to  reveal  the  array  frame. 
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Figure  8:  Airfoil  surface  pressure  (C;,),  Rcc  =  6.5  x  10^,  hjc  =  0.0019.  — o — *  —  0  (n  = 

0.020),  -  -X-  -  =  0.015  (ci  =  0.58),  --a--  =  0.057  {ci  =  1.5).  The  left  vertical  dashed 

line  represents  the  boundary  layer  trip,  and  the  right  vertical  dashed  line  corresponds  to  the 
slot. 


Figure  9:  Lift  coefficient  (q)  as  a  function  of  momentum  coefficient.  Rcc  —  6.5  x  10'^,  h/c 
=  0.0019.  o  present  investigation,  x  Abramson  (1975)  20%  elliptic  circulation  control  airfoil, 
Rec  =  3.4  X  10^  h/c  =  0.0013. 
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Figure  10:  Photograph  of  closed  test  section. 


Figure  11:  Airfoil  surface  pressure  (Cp),  =  0,  hje,  =  0.0019.  — o —  open  test  section, 

Re^  =  6.5  X  lO'*^;  -  -x-  -  closed  test  section,  Rvq.c  =  6.7  x  10^.  The  left  vertical  dashed  line 
represents  the  boundary  layer  trip,  and  the  right  vertical  dashed  line  corresponds  to  the  slot. 
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Figure  12:  Airfoil  surface  pressure  (Cp),  h/c  =  0.0019.  — —  open  test  section,  Cp  ~  0.057, 
Rcc  =  6.5  X  10^  (c/  =  1.5);  -  -x-  -  closed  test  section,  Cp,r  =  0.057,  Rccx  =  6.7  x  10^  (c/  =  2.7). 
The  left  vertical  dashed  line  represents  the  boundary  layer  trip,  and  the  right  vertical  dashed 
line  corresponds  to  the  slot. 


Figure  13:  Lift  coefficient  (q)  as  a  function  of  inomeiituin  coefficient,  Rr.c  =  6-5  x  lO'*^,  h/c  = 
0.0019.  o  present  investigation,  open  test  section;  A  present  investigation,  closed  test  st^ction; 
X  Abramson  (1975)  20%  elliptic  circulation  control  airfoil,  Rcc  =  3.4  x  10'^,  h/c  =  0.0013. 
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Figure  14:  Photograph  of  closed  test  section  for  PIV  measurements.  The  foam  ceiling  and  one 
foam  wall  are  replaced  with  clear  polycarbonate  panels. 


Figure  15:  Airfoil  surface  pressure  (C^)  in  closed  test  section,  Rcc  =  5.7  x  10^,  h/c  =  0.0019. 
— o —  foam  walls,  upper  slot  blowing,  =  0.057  (Q  =  2.7);  -  -x-  -  foam/polycarbonate  walls, 
upper  slot  blowing,  =  0.052  (Q  =3.1):  —  A--  foam/polycarbonate  walls,  lower  slot  blowing. 

=  0.053  {Cl  =  2.8).  The  left  vertical  davShed  line  represents  the  boundary  layer  trip,  and 
the  right  vertical  dashed  line  corresponds  to  the  slot. 
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Figure  16:  Airfoil  surface  pressure  (Cp)  in  open  test  section.  Rec  =  6.5  x  h/c  0.0019. 
— o —  upper  slot  blowing,  Cp  =  0.058  (c/  =  1.5);  -  -x-  -  lower  slot  blowing,  =  0.059  (c/ 
—  1.7).  The  left  vertical  dashed  line  represents  the  boundary  layer  trip,  and  the  right  vertical 
dashed  line  corresponds  to  the  slot. 


Figure  17:  Potential  flow  Cp  on  a  20%  ellipse,  U^c  =  20  in/s.  —  freestreani  (q  =  1.90); - 

-  ceiling  plane  (q  =  2.29); - ground  plane  (q  =  1.81);  o  o  c  ceiling  and  ground  plain's  (q  = 

2.26). 
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Figure  18:  Airfoil  surface  pressure  (Cp)  in  closed  test  section  compared  with  potential  flow 
theory,  x  foam/ polycarbonate  walls,  upper  slot  blowing  {Cci  =  3.1);  A  foam/polyearbonate 

walls,  lower  slot  blowing  (r/  =2.8); - potential  How,  ceiling  plane  (c/  =  2.29); - potential 

How,  ground  plane  (c/  =  1.81). 
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Figure  19:  Leading  edge  mean  flow  streamlines  in  (a)-(c)  open  test  S(x:tioii,  Rcc  =  G.5  x 
and  (d)-(f)  closed  test  section,  Rec,c  =  6-7  x  10^.  Flow  is  from  left  to  right,  and  the  origin  is 
located  at  the  leading  edge.  Lower  slot  blowing  is  used  for  all  cases,  /i/c  =  0.0019.  a)  =  0. 
b)  =  0.014.  c)  =  0.057.  d)  =  0.  e)  C^  c  =  0.013.  f)  =  0.053. 
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Figuro  20:  Comparison  of  profiles  at  the  lip  edge,  =  0.  •  open  test  section,  Rtc  =  G.5  x  10  ^ 
X  closed  test  section,  Rcc^c  =  6,7  x  10^.  The  slot  height  is  hfc  =  0.0019  for  both  cases.  For 
clarity,  only  every  fourth  datum  point  is  displayed.  Error  bounds  are  represented  by  lines  of 
matching  color  for  each  data  set. 
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F'igiiro  21:  Comparison  of  profiles  at  the  lip  odgo.  •  open  test  section,  Rtc  =  6.5  x  10^. 

=  0.014;  X  closed  test  section,  Rec^c  =  6.7  x  10^\  —  0.013.  The  slot  height  is  h/c  =  0.0019 

for  both  cases.  For  clarity,  only  every  fourth  datum  point  is  displayed.  Error  bounds  are 
represented  by  lines  of  matching  color  for  each  data  set. 
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Hgiiro  22:  Comparison  of  profiles  at  the  lip  edge,  •open  test  section  Re^  =  6.5  x  10'\  C^^  — 
0.057;  X  closed  test  section,  Recx  =  6.7  x  10^,  =  0.053.  The  slot  height  is  h/c  =  0.0019 

for  both  cases.  For  clarity,  only  every  fourth  datum  point  is  displayed.  Error  bounds  are 
represented  by  lines  of  matching  color  for  each  data  set. 
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Figure  23:  Logarithmic  spiral  surface. 


Figure  24:  Lengtli  and  velocity  scales  suggested  by  Launder  Sz  Rodi  (1983). 


Figure  25:  Coanda  surface  Cp.  Please  refer  to  Table  3  for  test  conditions  for  each  case. 


' 
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Figure  26:  Circulation  control  curved  wall  jet  profiles,  Rec  =  6.5  x  10'^,  hjc.  =  0.0019.  c^/c 
(a)  0.0036.  (b)  0.011,  (c)  0.018,  (d)  0.026,  (e)  0.033  (f)  0.040,  and  (g)  0.048.  •  C,,  =  0.015. 
=  0.0039,  C^,  =  0.057. 
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Figure  27:  Length  and  velocity  scales  based  on  zero  shear. 
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Figure  28:  Circulation  control  curved  wall  jet  profiles,  Rcc  =  6.5  x  10''’,  C,,  =  0.015.  h/<:  = 
0.0019.  rf)/c  =  (a)  0.0027,  (b)  0.010,  (c)  0.017,  (d)  0.025,  (e)  0.032  (f)  0.040,  and  (g)  0.017. 
Length  scales  are  indicated  by  o  ymhi,  x  Ve,  ^  o  ye,i/2- 
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Figure  29:  Outer  region  similarity  of  U  using  (a)  mean  shear  scales  and  (b)  scales  suggested  by 
Launder  k  Rodi  (1983),  Rcc  =  6.5  x  10^,  =  0.015,  h/c  =  0.0019.  rO/c  =  0.010,  a  0.017,  x 

0.025,  *  0.032  4*  0.010,  and  o  0.047.  The  line  represents  C/“0.5[(l  —  1.315  tanh(y“  1)]. 
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Figure  30:  Lack  of  inner  region  similarity  of  (7,  Rtc  =  6.5  x  10*'^,  =  0.015,  h/c  =  0.0019. 

rBjc  =  0.010,  A  0.017,  x  0.025,  ♦  0.032  -f  0.040,  and  o  0.047. 
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Figure  31:  Velocity  and  length  scales  describing  decay  and  spread  of  jet,  Rcc  =  6.5  x 
C/,  =  0.015,  and  h/c  =  0.0019. 


Figure  32:  Rate  of  spread  of  ymax/f^  and  collapsed  data  with  best  fit  line,  =  0.95.  Sec  Table 
3  for  test  conditions:  o  case  1,  case  2,  a  case  3,  x  case  4,  *  case  5,  o  case  6. 
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Figure  33:  Rate  of  spread  of  2/m,i/2/^^  collapsed  data  with  best  fit  line,  =  0.97.  See 
Table  3  for  test  conditions:  o  case  1,  c^ise  2,  a  case  3,  x  case  4,  *  case  5,  o  case  6. 


Figure  34:  Rate  of  decay  of  Umax/Ujet  and  collapsed  data  with  best  fit  line,  =  0.95.  See 
Table  3  for  test  conditions:  o  case  1,  case  2,  a  case  3,  x  case  4,  ♦  case  5,  o  case  6. 
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Figure  35:  Rate  of  decay  of  Uunn/Ujet  for  all  cases  listed  in  table  3,  and  collapsed  data  with 
best  fit  line,  =  0.95.  See  Table  3  for  test  conditions:  o  case  1,  case  2,  a  case  3,  >  case  4.  * 
case  5,  o  case  6. 


Figure  36:  Comparison  of  measured  and  predicted  (a)  velocity  scales  and  (b)  length  scales. 
Rcc  =  6.5  X  10'^,  Cfi  =  0.015,  h/c  —  0.0019.  Measured  scales  are  represented  by  datum  points, 
and  predicted  scales  are  represented  by  lines. 
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Figure  37:  Comparison  of  measured  and  predicted  (a)  velocity  scales  and  (b)  length  scales, 
Re-c  =  1.3  X  10^,  =  0.014,  h/c  =  0.0029.  Measured  scales  are  represented  by  datum  points, 

and  predicted  scales  are  represented  by  lines. 


Figure  38:  Mean  speed  (|V'|/f/oo)  contours  and  mean  velocity  vectors,  Re^  =  6.5  x  10'*^  and  h/c 
=  0.0019.  The  airfoil  trailing  edge  is  sketched  for  reference,  a)  =  0.  b)  Upper  slot  blowing. 
=  0.014.  c)  Upper  slot  blowing,  =  0.057. 
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Figure  39:  Approximate  separation  location  as  a  function  of  C^iRcc^  o  measured,  -  -  -  fit 
0.996). 


Figure  10:  Instability  range,  indicated  by  Gortler  number,  Re^  =  6.5  x  10'\  =  0.015. 

h/c  =  0.0019.  The  neutral  stability  curve  is  computed  by  Floryan  (1986). 


z/h 


Figure  41:  Instantaneous  spanwise  vorticity  (iOxh/Ujet)  distribution.  Re.c  =  5.6  x  lO'"’,  — 

0.014,  h/c  =  0.0019.  The  origin  of  z/h  is  located  at  midspan. 
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Figure  42:  Instantaneous  spamvise  vorticity  distribution  and  instantaneous  velocity 

vectors  reveal  presence  of  counter-rotating  pair  of  vortices,  Rcc  —  5.6  x  10^,  =  0.014, 

hjc  =  0.0019.  The  origin  of  z/h  is  located  at  inidspan. 
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Figure  13:  Power  spectral  density  measured  by  microphone  Ml  above  the  trailing  edge,  Re^  — 
6.5  X  10^,  h/c  =  0.0019. 
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-  /i/r  =  0.0010.  r  =  6.5  X  10^  =  0 

— - —  /i/r  =  0.0010.  =  I  .*!  X  10«.  =  0 

- hfc  =  0.00 10.  =  1.3  X  10®.  =  5  X  10- 

hfc  -  0.0010.  Rfr  =■  1.3  X  10®.  =  O.OOl 

- hfc  =  0.0019.  Rir  =  6.5  X  10\  =  0 

- -  h/c  =  0.0019.  Rtr  =  1.3  X  10®.  =  0 

- h/c  =  0  0019.  Re,  -  1.3  X  10®.  =  0.0111 

h/c  =  0.0029.  Re,  =  6.5  x  10^  =  0 

h/c  =  0.0029.  Rr,  =  13  x  10®.  =  0 

h/c  =  0.0029.  Re,  =  1.3  x  10®.  =  0.001 


Figure  44:  All  cases  with  trailing  edge  vortex  shedding  tones  measured  by  microphone  Ml.  The 
vertical  line  corresponds  to  St  =  0.21. 
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Figure  45:  Slot  lip  vortex  shedding  tones  measured  by  microphone  Ml,  =  0.014,  fijc  = 
0.0029. 
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Figure  46:  Slot  lip  vortex  shedding  tones  measured  by  microphones  Ml  and  M4  (see  Figure  4). 
Rcr  =  1.3  X  10^,  =  0.014,  h/c  =  0.0029. 


Figure  47:  Spectra  (16  Hz  bin  width)  measured  by  the  center  microphone  above  the  trailing 
edge,  SSB,  U^c  =  0,  fi/r  =  0.0019.  The  vertical  lines  represent  the  frequencies  where  beam 
maps  are  produced. 
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x/c 

Figure  48:  Sample  beam  map  with  airfoil  and  tunnel  components  labeled. 


Figure  49;  Outer  array  beam  maps  (dB)  at  992  Hz,  SSB,  Uo^  —  0,  fi/c  =  0.0019.  a)  Array  PSF. 
b)  Rcj.t  =  1340.  c)  =  2680.  d)  Rcjet  =  3980.  e)  Rvjet  =  5300.  f)  Rr^^t  =  6600. 
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Figure  50:  Outer  array  beam  maps  (dB)  at  2  kHz,  SSB,  Uoo  =  0,  /i/c  =  0.0019.  a)  Array  PSF. 
b)  =  1340.  c)  Re  jet  =  2680.  d)  Re^ei  =  3980.  c)  Rej^t  =  5300.  f)  Rejet  =  6600. 
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Figure  51:  Outer  array  beam  maps  (clB)  at  4  kHz.  SSB,  =  0.  h/c  =  0.0019.  a)  Array  PSF. 
h)  =  1340.  c)  Rej^t  =  2680.  d)  Rejf,t  =  3980.  e)  Rejpt  =  5300.  f)  Rejct  =  6600. 
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Figure  52:  Outer  array  beam  maps  (dB)  at  8  kHz,  SSB,  Uoo  =  0,  h/c  =  0.0019.  a)  Array  PSF. 
b)  Rejf>t  =  1340.  c)  Rejet  =  2680.  cl)  Rejet  =  3980.  e)  Rejet  =  5300.  f)  Rej^t  =  6600. 
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Figure  53:  Inner  array  beam  maps  (dB)  at  16  kHz,  SSB,  h/c  =  0.0019.  a)  Array  PSF. 

b)  Rej^t  =  1340.  c)  Rcjet  =  2680.  d)  Rej^t  =  3980.  e)  Rcjet  =  5300.  f)  Rejet  =  6600. 
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Figure  54:  Inner  array  beam  maps  (clB)  at  32  kHz,  SSB,  Ux  =  0.  k/c  =  0.0019.  a)  Array  PSF. 
b)  Rejrt  =  1310.  c)  Rejet  =  2680.  d)  Rejn  =  3980.  e)  Rej^,  =  5300.  f)  Re.jei  =  6600. 
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Figure  55:  Inner  array  beam  maps  (ciB)  at  64  kHz,  SSB,  Ux  =  0,  h/c  =  0.0019.  a)  Array  PSP', 
b)  Rcjrt  =  1340.  c)  Rcj^t  =  2680.  d)  Rejet  =  3980.  e)  Rej^t  =  5300.  f)  Re.j^t  =  6600. 
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Figure  56:  Spectra  (16  Hz  bin  width)  measured  by  the  center  microphone  above  the  trailing 
edge,  SSB,  Rcc  —  6.5  x  10^,  h/c  =  0.0019.  The  vertical  lines  represent  the  frequencies  where 
beam  maps  are  produced. 
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Figure  57:  Outer  array  beam  maps  (dB)  at  992  Hz,  SSB,  Rtc  =  G.5  x  10^.  h/c  =  0.0019.  a) 
Array  PSF.  b)  C,,  =  0.  c)  C„  =  0.004  =  1330).  d)  C,,  =  0.017  [Re^et  =  2660).  e) 

C,,  =  0.037  {Rcjet  =  4000).  f)  C,,  =  0.065  (Re^^t  =  5300).  g)  C„  =  0.10  (Re^^  =  6640). 
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Figure  58:  Outer  array  beam  maps  (dB)  at  2  kHz,  SSB,  Rep  =  6.5  x  lO’’,  h/c  =  0.0019.  a)  Array 
PSF.  b)  =  0.  c)  C„  =  0.004  =  1330).  d)  C„  =  0.017  {Re jet  =  2660).  e)  C„  -  0.0.37 

{Rejei  =  4000).  f)  Cf,  =  0.065  {Rcjei  =  5300).  g)  C,,  =  0.10  {Reje,  =  6640). 
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Figure  59:  Outer  array  beam  maps  (dB)  at  4  kHz,  SSB,  Rcr  =  6.5  x  10^.  h/c  =  0.0019.  a)  Array 
PSF.  b)  C,,  =  0.  c)  C,,  =  0.001  {Rejet  =  1330).  d)  C,.  =  0.017  {Rejei  =  2660).  e)  C,,  =  0.037 
{Rejei  =  4000).  f)  C„  =  0.065  {Rejet  =  5300).  g)  C,,  =  0.10  {Rcje,  =  6640). 
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Figure  GO:  Outer  array  beam  maps  (dB)  at  8  kHz,  SSB,  Rcc  =  6.5  x  10^,  h/c  =  0.0019.  a)  Array 
PSF.  b)  C,,  =  0.  c)  C,,  =  0.004  {Rej,t  =  1330).  d)  C,,  =  0.017  {Rejn  =  2660).  e)  =  0.037 
{Rej,,t  =  4000).  f)  C,,  =  0.065  =  5300).  g)  C,,  =  0.10  {Rcjet  -  6640). 
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Figure  61:  Inner  array  beam  maps  (dB)  at  16  kHz,  SSB,  Rcr  =  6.5  x  10’^,  h/c  =  0.0019.  a) 
Array  PSF.  b)  C,,  =  0.  e)  =  0.004  {Rejet  =  1330).  d)  C,,  ^  0.017  {Re^et  =  2660).  e) 
C„  =  0,037  {Rejet  =  4000).  f)  C„  =  0.065  {Rejei  =  5300).  g)  C,,  =  0.10  {Rejet  =  6640). 
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Figure  62:  Inner  array  beam  maps  (dB)  at  32  kHz,  SSB,  Re.c  =  6.5  x  10^,  h/c  =  0.0019.  a) 
Array  PSF.  b)  C„  =  0.  c)  C„  =  0.004  {Re^ct  =  1330).  d)  =  0.017  {Rej,t  =  2660).  e) 
C,,  =  0.037  {Rcj^i  =  4000).  f)  C,,  =  0.065  {Re^^t  =  5300).  g)  C,,  =  0.10  {Re jet  =  6640). 
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Figure  63:  Inner  array  beam  maps  (dB)  at  64  kHz.  SSB,  Re.^  =  6.5  x  lO'*,  hfe  —  0.0019.  a) 
Array  PSF.  b)  C^,  =  0.  c)  C,.  =  0.004  =  1330).  d)  Q,  =  0.017  {Re^^t  =  2660).  e) 

=  0.037  {Rej^t  =  4000).  f)  C,,  =  0.065  [Rcjei  =  5300).  g)  C,,  =  0.10  {Re.jrt  =  6640). 
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Figure  61:  Outer  array  beam  maps  (clB)  at  992  Hz,  Re,.  =  6.5  x  10^,  h/c  =  0.0019.  a)  SSB. 
C„  =  0.037.  b)  OSB,  C,,  =  0.040.  c)  SSB,  C„  =  0.065.  d)  OSB,  =  0.071.  e)  SSB,  C„  =  0,10. 
f)  OSB,  C„  =  0.11. 
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Figure  65:  Outer  array  beam  maps  (dB)  at  2  kHz,  Rcc  =  6.5  x  10'’,  h/c  =  0.0019.  a)  SSB. 

=  0.037.  b)  OSB,  C,,  =  0.040.  c)  SSB,  Q  =  0.065.  d)  OSB,  C„  =  0.071.  e)  SSB,  =  0.10. 
f)  OSB,  Q  =  0.11. 
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Figure  66:  Outer  array  beam  maps  (clB)  at  4  kHz,  Rec  =  6.5  x  lO'^,  h/c  =  0.0019.  a)  SSB, 
C,,  =  0.037.  b)  OSB,  C,,  ^  0.040.  c)  SSB,  C,,  =  0.065.  d)  OSB,  C,,  =  0.071.  e)  SSB.  C„  =  0.10. 
f)  OSB,  C,,  =  0.11. 
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['igure  67;  Outer  array  beam  maps  (clB)  at  8  kHz,  Rcr  =  6.5  x  lo'*,  h/c  =  0.0019.  a)  SSB. 
C,,  =  0.037.  b)  SSB,  C,,  =  0.040.  c)  SSB,  C,,  =  0.065.  d)  OSB,  C,,  =  0.071.  e)  SSB.  =  0.10. 
f)  OSB,  C„  =  0.11. 
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Figure  68:  Inner  array  beam  maps  (dB)  at  16  kHz,  Rcc  =  6.5  x  10^,  li/c  =  0.0019.  a)  SSB. 


=  0.037.  b)  OSB,  =  0.040.  c)  SSB,  =  0.065.  d)  OSB,  C,,  =  0.071.  e)  SSB,  =  0.10. 
f)  OSB,  =  0.11. 
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Figure  69:  Inner  array  beam  maps  (dB)  at  32  kHz,  Re.c  =  6.5  x  10',  h/c  =  0.0019.  a)  SSB, 
C„  =  0.037.  b)  OSB.  ==  0.040.  c)  SSB,  T,,  =  0.065.  d)  OSB,  C,,  =  0.071.  e)  SSB.  T,,  =  0.10. 
f)  OSB,  =  0.11. 
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Figure  70:  Inner  array  beam  maps  (dB)  at  04  kHz,  Re^  =  6.5  x  10^,  h/c  =  0.0019.  a)  SSB, 
C„  =  0.037.  b)  OSB.  C,,  =  0.040.  c)  SSB,  C,,  =  0.0G5.  d)  OSB,  C„  =  0.071.  e)  SSB,  Q,  =  0.10. 
f)  OSB,  C,.  =  0.11. 
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Figure  71:  Compcirison  of  free-standing  iiiicrophoiie  processing  techniques  for  SSB,  = 
6.5  X  10'^,  =0.11,  and  h/c  =  0.0019  (16  Hz  bin  width). 


Figure  72:  Comparison  of  free-standing  microphone  processing  UKihiiiques  for  OSB,  =  0.10, 
Re^  =  6.5  X  10^,  and  h/c  =  0.0019  (16  Hz  bin  width). 
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Figure  73:  Single  microphone  and  integrated  array  spectra  (16  Hz  bin  width)  for  SSB,  % 
0.065,  Rec  =  6.5  x  10^  and  h/c  =  0.0019. 
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Figure  74:  Mean  tangential  velocity  profile  along  trailing  edge  at  rO/c  =  0.0129,  Rec  =  6.5  x  10^"*, 
=  0.057,  h/c  =  0.0019  with  one-seventh  power  curve  inset  (Praadtl  1961). 
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Figure  75:  Meaui  tangential  velocity  profile  at  slot  lip  edge,  Rcc  =  6.5  x  10^,  =  0.057,  fi/c 

=  0.0019.  Notice  the  local  maximum  velocity.  Only  every  other  point  is  displayed  for  clarity. 


Figure  76:  Mean  tangential  velocity  profile  at  slot  exit  with  one-seventh  power  curve.  /?c,.  = 
6.5  X  10\  =  0.057,  h/c  =  0.0019  (Prandtl  1961). 
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Figure  77:  Howe's  model,  evaluated  for  Rcr  =  6.5x  10‘\  =  0.057,  h/c  =  0.0019,  providing  only 

airfoil  geometry  details  and  test  conditions.  The  array  spectrum  represents  SSB,  Rcc  =  0.5  x  10^, 
Cf,  =  0.065,  h/c  =  0.0019. 
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Figure  78:  Howe’s  model,  evaluated  for  Rcc  =  6.5  x  10^\  =  0.057,  h/c  —  0.0019  using 

length  and  velocity  scales  estimated  from  PIV  data.  The  array  spectrum  represents  SSB. 
Rcr  =  6.5  X  10^  C,,  =  0.065,  h/c  =  0.0019. 
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Figure  79:  Best  fit  of  Howe's  model  to  measured  speetruiii. 


Table  1:  Scaling  of  length  and  velocity  scales  from  similarity  solution  of  Guittoii  k  Newman 
(1977);  is  the  arc  length  and  n  is  a  constant 


Parameter 

Scaling 

R 

X 

.V 1/2771 

X 

TJ 

TTiaO* 

Table  2;  Scaling  of  length  and  velocity  scales  from  general  similarity  solution  of  a  curvenl  wall 
jet  in  a  fret^stream;  x  is  the  arc  length  and  a  is  a  constant 


Parameter 

Scaling 

R 

X 

!J\ 

./■ 

\J2 

Ih 

j" 

U2 

x" 

Ih/lh 

constant 
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Table  4:  Length  and  velocity  scales  usee 
Predicted  (Howe  2002) 

U  (ni/s)  (mm)  u*  (m/s) 


Noise  type 


Curvature 
Passive  slot 
Slot-jet 


56.3 

32.0 

79.9 


1.04 

1.04 

0.019 


1.97 

1,12 

2.80 


to  evaluate  Howe's  model 
PIV  Measurement 
U  (m/s)  (mm)  lu  (m/s) 


91.3 

41.4 

88.8 


0.010 

0.140 

0,025 


9.57 

2.78 

6,84 


Table  5:  Example  of  scales  required  for  Howe's  model  to  match  measurement 


Noise  type 

U  (m/s) 

6^  (mm) 

v,  (rn/s) 

Curvature 

63.1 

0.5 

20.0 

Passive  slot 

59.8 

0.17 

19.0 

Slot-jet 

88.8 

1.0x10-^ 

7.60 
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